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PART A. INCOMPRESSIBLE PLOW 

1. Introduction 

The purpose of the stability theory of the laminar boundary layer is 
to determine whether a small disturbance introduced into the boundary layer 
will amplify or damp. If the disturbance damps, the boundary layer should 
remain laminar. If the disturbance amplifies, and by a sufficient amount, 
then transition to turbulence should eventually take place. The stability 
theory cannot predict the location of transition. What it can do is to 
establish those states of the boundary layer which are most likely to lead 
to transition, - to identify those frequencies which are the most dangerous, 
and to indicate how the external parameters can best be changed to avoid 
transition. 

It’ will be the purpose of these lectures to give the fundamentals of 
the stability theory, to establish its chief results, and to give an idea 
of the physical mechanisms at work. Detail, except for a few topics felt 
to be of particular importance, will be left to the numerous references. 

Of the latter, there now exist a number of review articles and two books. 
The articles are by Shen (1964-), Stuart (1963), Reid (1965), and Drazin and 
Howard (1966). The latter is devoted exclusively to the inviscid theory. 
The article by Reid is particularly complete on the asymptotic theory. The 
books are by Lin (1955) and Betchov and Criminale (1967). An older review 
article that can still be read with profit is Prandtl (1935)* The well- 
known book of Schlichting (i960) on boundary-layer theory includes two ■ 
chapters on boundary-layer stability. 
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There has always been a certain air of mystery about stability theory 
for many people. This is probably a result of the elaborate mathematical 

procedures that were necessary to arrive at any results. The well-developed 

• 

inviscid stability theory suffered the embarrassment of predicting complete 
stability for two of the inoet common flows , the Blasius boundary layer and 
flow in a pipe. The viscous theory was of such a complexity that almost 
no general results could be established. In spite of the many difficulties, 
the instability of the incompressible Blasius boundary layer was worked out 
by Tollmien (1929), and this theory received confirmation in its essential 
aspects in the classic experiment of Schubauer and Skramstad (1947). 

Since about i960, the approach to stability theory has been considerably 
modified by the modern high-speed digital computer. Instead of the painful 
extraction of a handful of numerical results from elaborate and time-consuming 
calculations of dubious accuracy, the computer can produce an immense quantity 
of accurate numerical results directly from the differential equations. The 
ease of producing these results makes it possible to not only solve any specific 
problem, but also, to fill in the gaps in our understanding which still exist 
because of the difficulty of analysis. For example, the theory deals almost 
exclusively with disturbances of neutral stability, which from a practical 
standpoint are of limited interest. The computer permits the study of the 
much more interesting disturbances of maximum amplification, which can be 
obtained as easily, indeed more easily, than the neutral disturbances. 

In spite of the limitations of the theory in producing numerical results, . 
it is still essential for our understanding of boundary-layer stability and 
as a guide to the formulation and use of the computer programs. For this 
reason, it will be covered in these lectures, and will provide the motivation 
for the entire subject. The computer will elaborate and provide the main body 
of numerical results. 
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Before we get. into the main body of the subject, a few words are in 
order to orient those who are new to this field. The stability theory 

considers individual sinusoidal disturbances whose amplitude is small 

• * 

enough so that a linear theory can be used. The wave number in the free- 
stream direction is nr/\ > where X is the wave length). The 

wave may be two-dimensional, with the ' normal to the wave front parallel to 
the free-str'eam direction, or it may be oblique, with the wave normal at 
an angle \js to the free- stream direction. The disturbance propagates in • 
the downstream direction with phase velocity Cj. . The phase velocity is 
always less than the free-stream velocity, "Uj , so that at some point in 
the boundary layer the mean velocity is equal to c*. . This point is called 
the critical point, and it plays a central role in the mathematical theory. 

At a given distance from the origin of the boundary layer, or better, 
at a given Reynolds number R* ( - 1 1,%/ir 9 where u- is the kinematic 
viscosity), the disturbance may be in one of three states: damped, neutral, 

or amplified. The numerical results obtained from the stability theory are 


often presented in the form of neutral-stability 'diagrams which show graphi- 
cally which state a given disturbance is in at each Reynolds number. There 

t 

are two general kinds of neutral- stability diagrams to be found, as is shown 
in Fig. 1.1. 



Fig. 1.1, Typical curves of neutral stability 
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In these diagrams, is plotted against R* . Disturbances are neutral 
at those values' pf and which lie on the contour marked neutral; they 
are amplified inside of the contour; and they are damped - putside of the 
contour. With a neutral-stability curve of type (a), ‘all wave numbers 
are damped at sufficiently high Reynolds numbers. Xn this case;, the mean 
flow is said to have invlscid stability. Since decreasing Reynolds number, 
or increasing viscosity, can lead to instability, it is apparent that vis- 
cosity' does not act solely to damp out disturbances, but actually can have 
a destabilizing influence. The flat plate, or Blasius, boundary 'layer is 
an example of a flow which is unstable only through the action of viscosity. 
With a neutral-stability curve of type (b), wave numbers smaller than c<y > 
the neutral wave number at infinite Reynolds number, remain unstable no 
matter how large the Reynolds number becomes. A mean flow with a type (b) 
neutral- stability curve is said to have inviscid instability. The boundary 
layer in an adverse pressure gradient is an example of a flow of this kind. 

In both cases (a) and (b), all disturbances with ®{ < are unstable 

for some range of Reynolds numbers. The Reynolds number, Rei- below which 
no amplification at all is possible, is called the critical Reynolds number. 
It is often an objective of the stability theory to compute f?c,- , although * 
it must be cautioned that it has a limited significance and is often not the 
best indicator of the relative instability of various mean flows. 

Although the wave number is a useful quantity in developing the theory, 
the frequency, which is equal to cx c t . , is more useful in practice. A dis- 
turbance which is introduced into the boundary layer v-ith a particular fre- 
quency will preserve that frequency as It progresses downstream;- 'while the 
wave number will change. A frequency f which crosses the unstable region 
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will "be damped from the leading edge up to Rl , the first neutral point. 
Between and R v , the second neutral point, it will "be amplified; 
.downstream of f? v it will he damped again. If the magnitude of. the dis- 
turbance' becomes large enough before Rv is reached, then' the nonlinear 
processes which eventually lead to transition will take <?ver, and the 
disturbance will continue to grow even though the linear theory says it 
should damp. 

The theory can be used to calculate amplification and damping rates 
as well as the wave number, Reynolds number and phase velocity of neutral 
disturbances. For example, it is possible to compute the amplification rate 
as a function of frequency at a given R* . The neutral -stability curve only 
identifies the range of unstable frequencies, but the calculation of the 
amplification rates tells how fast each frequency is growing, and which 
frequency is growing the most. Perhaps even more useful than the ampli- 
fication rate is the overall growth of a disturbance of constant frequency 
as it travels through the unstable region. This result can be obtained from 
the theory as a supplemental calculation once the amplification rates are 
known. Consequently, it is possible to identify, for each initial spectrum, 
the frequency which has the largest amplitude at each Reynolds number. It _ 
is presumably one of these frequencies which, when it reaches some critical 
amplitude, triggers the whole transition process. - ' 
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Derivation of Stability Equations 


2.1 Formulation of theory 

i 

The conventional form of the theory will be used throughout these 

* • 

lectures. That is, the actual boundary layer is replaced by, a- parallel 
flow having the same velocity profile as the boundary layer at a particular 
Reynolds number. A single Fourier component is considered as the disturbance, 
Xt extends to Infinity both upstream and downstream and grows or damps in 
time, not in space. However, the physical problem suggests a disturbance 
growing in space, and the theory could he carried through for this type 
of disturbance. 


The 

Navier* 

•Stokes equations 

of a viscous 

incompressible fluid are: 
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The x axis can be considered to be in the direction of the free-stream 
velocity j the y axis is normal to the surface; and the z axis is normal 
to the x and y axes. The velocities u , \r , w are in the x, y, z 
directions, respectively. The density is p , the pressure p , and the 
kinematic viscosity a . The first three equations are the x, y, 
z momentum equations; the fourth equation the equation, of continuity. 
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2.2 Linearization of equations 

All flow quantities are divided into a mean flow ter^i and a fluctuation 


term: 
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(£. 5 ) 


The mean-flow terms (denoted by capital letters) satisfy the boundary-layer 
equations^, which are 
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In the linear case it will not be necessary to distinguish between 
the undisturbed mean flow and the mean flow in the presence of a disturbance. 
The replacement for (2.6) in the nonlinear case is given in Section 10 by ■ 
(10.6), When the expressions (2,5) are substituted into (2.l)-(2.4); the 
negligible mean -flow terms (by the boundary-layer approximation) dropped j 
the mean boundary-layer equations subtracted out$ and the nonlinear terms 
neglected j there remains the following system of equations: 
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2.3 Reduction to parallel -flow equations 

The above equations are still too complicated and must be reduced 
further. If it is assumed that all velocity fluctuations are of the same 
order of magnitude, and all derivatives of velocity fluctuations are of 
the same order of magnitude (but possibly different from the fluctuations), 
then the application of the boundary-layer relations to the mean flow 
quantities will result in the desired simplification. These relations 
are ? 
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The terms with velocity fluctuations are compared separately, and those 
with derivatives of the fluctuations are also compared separately. The 
equations reduce to 
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Now a close inspection reveals that the coefficients of the V terns in 

(2.15) and (2.I7) are an order of magnitude larger than the coefficients 

of the and 3/32: terms on the right-hand sides,, while in (2.16) 

the coefficient of V* is of the same order as the other two terms. Since 

\r , <.£ ILV for the linearisation to be valid, the v # 3Y/3 * 4 term can be 
* . 0 

dropped from (2.16). In (2.15) and (5,17), although 3 U/<?^ and OYr/Si 

and t - ' 9^/3^ 

terms are not necessarily larger than the terms immediately preceding and 
following, because although \r' is always small, the derivatives of lA * and 
V* do not have to be small. Consequently, (2,16) is replaced by 

!?' 4 ld'i -w L v ' = - 1 3 J' + * % 

3* 3iC H l 3^ (2.16a) 


are an 'Order of magnitude larger than U. and Yf , the v‘ 3 ZU 3 


and (2.15) and (2.17) remain as given. The final approximation is to con- 

1 

sider that for disturbances for which the wavelength is not too long, the 
mean-flow quantities will be slowly \ r arying functions of x and £ compared 
with the fluctuations. Therefore, V. and W can be taken as functions of 
o' alone, and the stability of the velocity profile at a. given x is con- 
sidered independently of the rest of the boundary layer. 
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Ihe final system of equations (2.15 ), (2.16a); (2.17); (2.18) is the 

same as would have been obtained by malting the. as sumption of parallel flow 

/ 

at the outset. This assumption is 


1C n llifi , 'WT - *WC^| , V = o 


(2,19) 


The argument given above that led to these equations is far from rigorous 
and is included mainly to point out the terms that are neglected in the 
parallel-flow equations. These equations are the basis of almost all 
stability investigations. They are exact for the flow in a channel, but 
are only an approximation for other flows. Their adequacy is best tested 
by comparison of the results with experiment. Also the magnitudes of the 
neglected terms can be checked a posteriori, and it is hoped that some day 
a more exact theory will be available for . comparison, 

2.4 Jfondimensional equations 

Before proceeding further, the equations will be written in dimen- 
sionless form. The characteristic velocity is Id, , the free -at learn velocity 
in the X direction; the characteristic length is S , the boundary-layer 
tkicknes’s; and the characteristic pressure is ip* , the free^stream pressure. 
If the dimensional quantities are designated by asterisks (except for the' 
boundary-layer thickness), the nondimen sional quantities are 
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and the Reynolds number is defined as 



( 2 . 21 ) 


We define the boundary-layer thickness to be equal to at the point where 
V? ~ 0.999 Uj . For the Blasius boundary layer Slif* / K* 53 6.0 t The 
dimensionless equations are the same as (2.15)-(2.l8) except for the coef- 
ficient of the pressure terms and the replacement of Jr by */[?$ . The 
pressure is nondimens ionalized with respect to the free -stream pressure 
rather than the more usual in order to make the development here 

correspond more closely to the derivation for compressible flow to be given 
later. * The dimensionless equations are 
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2,5 Boundary conditions 

The boundary conditions are that the no-slip condition applies to the 
disturbance velocities at the wall, 




0 f Co } - Q f w'Coj - 0 


(2.26) 
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and that the disturbances go to zero as ijVcj , 

-v u / ~v o t -? 0 as ^ -9 «*? (2.27) 

Since all of the boundary conditions are homogeneous, it can be expected 
that solutions will exist only for particular combinations of 1?^ and the 
parameters of the disturbance ( etj, , , c 3 see next paragraph). That 

Is, the stability problem is an eigenvalue problem. The values of otj j • 
^ ) c and for which the boundary conditions can be satisfied are 
eigenvalues,, and the corresponding amplitude functions are eigenfunctions. 
2.6 Introduction of Fourier components 
. The final form of the differential equations , where the coefficients 
are functions only of ^ , and * , £ 3 rb appear only as derivatives, 

suggests the following ^type of disturbance. 1 ... 

e-sp[iC«sX+ ( (2.28) 

In (2.20), is a typical complex amplitude function; and are the 

1 ' 

dimensionless wave numbers, - ZtT 5 f \\ , ^ - ItiS/V*. , where X* and 

1 j 

* * 

A 4 are the wavelengths in the >< and £• directions, respectively, and 
S is the boundary-layer thickness. If and {3 5 are real, then a 
complex c will give a disturbance that grows in time. With 

C ~ + id; • (2.29) 

the phase velocity is c w . The dimensionless frequency Oj is w ic : r ., and 
e*„ c* is the time derivative of the logarithm of the amplitude. 
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i- dll! = «t t c, (2.30) 

We shall refer to ot s C{ as the time rate of amplification. It is also 
possible to take c( s complex and c real* thus providing' a disturbance 
that grows in space. All results to be presented here will be for real 
<% and complex c » In (2.28)* only the real part of <j is to be con- 
sidered' to have physical meaning. 

Realij a. C ’ i cos -s- 5 - 0 < s Cf } 1 ( 2 . 31} 

5i h ( «5 X 1 Cir £)] 

Equation (2.30) gives the time rate of amplification* but we also 
need to know the space rate of amplification in order to compute the 
overall growth of a disturbance as it travels thrbugh the boundary layer. 

The space rate can be obtained from the time rate’ only if the proper prop- 
agation velocity is known. When a wave packet made up of individual waves 
of the form 

escf> i C <* x - wi) 

propagates through a dispersive medium* i.e.* a medium* such as the boundary 
layer* where the frequency Co is a function of cs( * each frequency advances 
with the phase velocity ' " 

c <r - 7 (2.32) 
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but an overall quantity such as the energy of the wave packet advances 
with the group velocity 
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It seems reasonable that the propagation velocity needed to obtain a space 
rate of amplification from (2.30) would be the group velocity rather than 
the phase velocity. Therefore^ we write 
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and the ratio of the amplitudes at two Reynolds numbers is 
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vft, c 3 

where the integration will normally be carried out for a constant dimen- 
sionless frequency, A commonly used form of the dimensionless frequency 
which varies as the frequency for a constant free stream is • 
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Equation (2.34) has teen demonstrated by Gaster (1963) to be correct 
for small amplification rates. Unfortunately the group velocity requires 
some additional calculations to be made* and it is difficult to resist the 
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temptation to simply replace it with the phase velocity. This temptation 
is all the more attractive in that C*. rarely differs from c** by more 
than about 1 5 $, and, perhaps more important, excellent agreement is 

obtained with experiment by using c v a Another difficulty with the group 

.s ' 

velocity is that if there is a region of anomalous dispersion ( c?^/ Ad 
changes sign), the group velocity no longer represents the velocity of 
energy propagation in this region* In the compressible stability theory, 
we will find that just such regions can actually occur* 
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3. Reduction of Equations to Two-Dimensional Form 

Inspection of (2.22)-(2.25) reveals that w* and "V ocpur only in 
9/32 terms except In the z momentum. equation* Consequently, if the 
disturbances have the form (2.28), and the x, z axes are rotated about 
. the y axis so that the new z axis is parallel to the wave, front, the z 
derivatives in the new coordinate system must be zero, and the 'and 
W terms will drop out of (2.22), (2.23) and (2.25). This special 
coordinate system is called the tilde system. If, further, the reference 
velocity is changed from to lC - Ut Cos ¥ , where, iff is the angle 
between the wave normal and the x axis, then the renormalized equations 
in the tilde coordinate system are 
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and 


U* * U' co s ¥ 




The original and tilde coordinate systems are shown in Fig* 3.1* 
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Fig, 3*1® Original and tilde coordinate systems 

. nta (&*£•<*£&> F? ^£> ? 

The system of equations (3*l)j (3*2) and (3**0 is)identical in 

respects to the equations for a two-dimensional' disturbance in the 

r+ /*jL A * 

boundary layer lUy\ at Beynolds numbe ^Ky, >-ghe transformation 
formulas are ' ■ 1 . ! 


U - 1L * “W -to* t 

%i - w* , w' 


~ t 


V' 1 / cos "f 


t _ 


“ £ COS ^ 


~Kr * - U **W v +■ " VJ " 

W ? - - u' *Kv\ 4- W f 
fc - - V 5n\\1f ■*» l CffS'vJ*' 


The dimensional normal velocity is. unchanged hy the transformation, 

- • i .. 

the dimensionless normal velocity must change because of the change 
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the reference velocity. The transformation is, cf course, only valid 
for disturbances of the type (2.28) with a wave front at a fixed, : angle. 
The relations between the wave numbers and wave velocities in the two 
coordinate systems are 

c< - C^ 2, t Icosty , 'c - c (3*8) 

i 

Although the dimensional phase velocity in the x direction is c* Cos \\/ , 
the dimensionless phase velocity 'c is the same as c because of the 
difference in the reference' velocities in the two cases, 

A* 

Since is different for each ty* , a different eigenvalue problem 

must be solved for each 'Vp* » The only simplification is that la lower 
order system of differential equations can be used. . However, to' obtain 
the disturbance velocities u 1 and V* in the original coordinate system, 
it is necessary to know w ^ , which means that the z momentum equation 


aw' ^ • 4 a± . i (isl : t t£! ) 

f 9* -l-J \ \ it »n‘7 


(3-9) 


■i: , .r 

must' be solved. Note that because of the term t \r , there will 

•v / ^ 

always be a w component as long as there is a."W* component. A normal 
fluctuation velocity always gives rise to a fluctuation in any^ quantity 
which has a normal gradient. * ( , 

If “W*- o in the original coordinate system, , i.e„, the 

A-* 7 

boundary layer is two-dimensional, 'll is identical to h , and t£e -brans- 

v 

formed equations for a three-dimensional disturbance at angle are 

, « 

precisely the equations for a two-dimensional disturbance at a lower 
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Reynolds number ce s , This is the content of the celebrated theorem 

of Squire (1935)* Consequently, if the minimum Reynolds number for in- 
stability of two-dimensional disturbances is Rev , then it follows that 
the minimum unstable Reynolds number for three-dimensional disturbances 
at wave angle is Rcv/dus^. It is in this sense that the statement 

is made' that the most unstable disturbance is a two-dimensional disturbance. 
Obviously, if the complete stability characteristics are known for two- 
dimensional disturbances, then the stability characteristics of all oblique 

waves can be obtained. For a wave at angle ^ and Reynolds number R 5 with 

wave ^ 

numbers o(<j and , the time rate of amplification is equal to ^ C\ cps\|/, 

where is the wave number of a two-dimensional disturbance at 

Reynolds number R s - . It is possible, depending on the exact 

shape of the unstable region in an ^ j R 5 plane \ that an which is 

stable for Mf-cO is unstable for some range of non-zero \Jr . Since the 

stability characteristics of a three-dimensional wave can be deduced from 

those of a two-dimensional wave, the subsequent analysis will be for a two- 

dimensional wave* 





4. Special Forms of the Stability Equations 
4.1 Orr-Sommerfeld equation 

The Fourier components of the flow quantities are defined to be 

-c4 
(4.1) 

f'Cx, y, t\ - TTCj) *» Cx-c;j6) 


u’u, ^ , i) = -fC^ s' 


« *6 (tf- C =£) 


V c* f ^ &) - ot £ 


i«; 


(if ip is replaced by , then ^ is the amplitude function of the 

stream function. This procedure is the usual one., but the present 
derivation is chosen to correspond with the course followed in com- 
pressible flow.) Substitution of these relations into (2.22), (2.23) 
and ( 2 . 25 ) yields 


IU-cH + U'j 


3 *rr + L 






(r- <\D 


$ - - -~7 1 ' ir ‘ + —(<$"- 'tjfi 
if + - o 


(4.2) 


(4.3) 


(4.4) 


The primes now refer to differentiation with respect to ^ . If (4.2) is 
differentiated with respect to ^ and 7T / and 4 eliminated, a single 4th- 
order differential equation is obtained for <p . 

Ck-Ofy"- < 1 $) - it y -- -- 1 — 0 ? lv - issf 1 - ^ ?) (4.5) 

This is the Orr-Sommerfeld equation. The boundary conditions are 
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J(o) - 0 

$ -v 0 


j'Ce) = 0 ■ 

o <vs Aj -V 


(4.6) 


t>6 


4. 2 Four first "Order equations 

it is also possible to write (4.2)-(4.4) as four first-order 
equations in terms of -f , -p' ; TT and ip . In order to correspond with 
the compressible -flow theory , these are written as 

_ K . 


*' , h - * ... * 


•4- - 




■t; 




The four equations are 


- z 


il - - i 2, 


l' -- - Li* i t - 4 
ft 


with the boundary conditions 


i 01- c) t ^*U 3 
. R« 


2, (o) ^ o t Zz Co) - <? 


-V ° j ^5^1 ’f’ 0 \j. 


(4.7) 


(4.8) 


i% * [* f? s ( 1t-c) + ^ ] 2i + <*5 $s & 2 3 + i oCtf f?s £* (4.9) 


(4.10) 


(4.11) 


(4.12) 
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4,3 Simplified forms of the Orr-Sommerfeld equation 

4.3.1 Fr.ee stream . 

The Orr-Sommerfeld equation will now be applied to three simplified 
situations: to the free -stream, where It - i , U! - 0 f 1 1 - 0 ; to 
inviscid disturbances, where s- w ; and to viscous layers within 

the boundary layer. 

In the free stream, (4*5) can be written 

1 




(4.13) 


Since this equation has constant coefficients, the solutions must have the 

form 


?<3) = Ae 


x 3 


(4.14) 


When this form is substituted into (4.13), a quartic equation for the 
characteristic value X is obtained (the term eigenvalue is reserved for ' 
those combinations of D<$ , (?$ , c which provide solutions of the Orr- 
Sommerfeld equation and boundary conditions). 


i « 4 I?, ({- c) (\ l - «s ) - (\ l - *» ) 


(4.15) 


The four solutions can be written down <•'/ inspection, and are 

X = 1 , X = i [ojf + i«sSiCl'tl] 1 


(4.16) 


Only two of these characteristic values satisfy the boundary conditions 
at infinity. They will be called and \i » 
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\ ~ .<*S ) X 3 - - [<i + iot s { 2 5 a-c^] (4.17) 

.* 

The corresponding characteristic vectors A are, from (4.8)^ (4.9)/ 

I 

(4.10)/ (4.1l)/ and leaving out A$^ / 

A W' 

A 1 = - 1 

. (4.18) 

A 5 - i 

A« = i«5 <l-c\ 

and 

( a \ r 2. i 

>. A, - - I [<*& + I *s Ri 

- l (4.19) 

A c>\ 

Aa ” 0 

, tt\ 

The vector A is the inviscid solution/ and is just the linearized inviscid 

, A <tt 

solution over a wavy vail of wavelength 2JT/c(i . The second vector A is 
the viscous solution, and represents a viscous decay. 

4.3.2 Inviscid equation 

The inviscid solution in the free stream is obtained from the inviscid 
equation 

Cj>“ ' <J - O ■ (4.20) 
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(4.21) 


which is the free-stream form of the complete inviscid equation 
— 

J ZC-h 7 

(4.21) is obtained from the complete Orr-Sommerfeld equation by taking .... 
the limit ^-* 00 . 

Since (4.21) is only a second -order equation and the Orr-Sommerfeld 
equation is fourth order, only two instead of four boundary condition® 
can be satisfied. These two conditions are 

- 0 7 o Aj -* (4.22) 

As is usual in inviscid flow, the normal velocity at * j - o is zero, but 
the no-slip condition is not’ satisfied. 

4.3.3 Simplified viscous equation . 

Since instability occurs at large Reynolds numbers, <x i will be 
large and, consequently, 

r -A 

X, a - [i« 4 Ks (4.23) 

It is interesting to observe that the version of the Orr-Sommerfeld 
equation that leads to (4.23) is 

(lt-c) <p" * - ' (4.24) 

T <*«(?» J 

Both (4.2l) and (4.24) are the basis for important further developments. 
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5- Inviscid Theory 


The inviscid equation, (4.21), is important for two reasons. 

First, it is possible that there are situations where viscosity is of 
little importance in the stability problem, except as it establishes 
the mean flow, and a purely inviscid stability theory is applicable. 
Second, the asymptotic viscous theory, where Rj is large, uses 
the inviscid solutions as two of the necessary four independent solutions. 
It is important to clearly distinguish between these two separate uses 
of the inviscid solutions. Since the inviscid equation is much simpler 
than the Orr-Soramerfeld equation, an extensive inviscid theory has been 
developed, mainly by Rayleigh (see list of references in Betchov and 
Criminale (1967) ), with important additions by Tollmien (1935) Lin 
(1945). A comprehensive review of the inviscid theory may be found in 
the article by Drasin and Howard (1966), 

5-1 Rayleigh theorems 

As a first step in the inviscid theory, some important results of a 
general nature will be established. Multiply (4.2l) by , the complex 
conjugate of <p . 

— H" . (1 x-tU <?<?* » 0 (5.1) 

Subtract from this equation its complex conjugate. Then 
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When (5*2) is integrated from n = 0 to infinity, the first term is 
zero by the boundary conditions. Hence 

c -, [ - -1 $ -- 4, = o ■ (5-3) 

Ifc-d > 

It follows from (5*3)’ that if c ; £.0 ,H must change sign in the interval 
0 < a t . Consequently., a necessary condition for an amplified dis- 
turbance is that the boundary-layer velocity profile must have an Inflection 
point. This result was first obtained by Rayleigh. It has subsequently 
been proved by Tollmen (1935 ) that for certain velocity profiles, which 
include the boundary-layer profile, this condition is also sufficient. A 
stronger result was obtained by Holland (1953)- He proved that an amplified 
solution can exist only if ZfY is negative between ^ 

Consequently, * the vorticity U' must have a maximum at * This condition 
rules out certain velocity profiles for which Ut-o , but where is a 
minim ?m„ 

It can also be proved that an inflection point is necessary for a 
neutral disturbance. Let 

= j'yM (5 ’ 4) 

It follows from (5*2, that with c ; =o,'W' must be constant except possibly 
at the critical point. A result of Rayleigh for which the proof will not 
be given establishes that the critical point will always lie between 
the wall and free stream. When (5*2) is integrated across the boundary 
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layer* the only contribution to the integral for C; -v o comes from 
near ^ • . 


\T(*j c + "VfC^j* - fil 


^ I?/ !'•« [ Alt 

K «-] (U-Cr) l f c! ■ 


( 5 . 5 ) 


The integration variable has been changed from aj to U . In the limit 
of q o , the integral of (5-5) acts like a delta function. Conse- 
quently, 


+ c) - e\ - IT -1 l^l 


( 5 . 6 ) 


Since by the boundary conditions, "WC^ c +£\ and are both zero, 

Vc must also be zero, and it has been" proven that a neutral disturbance 
can only exist if the velocity profile has a point of inflection. It 
further follows from (5*5) that the phase velocity of the neutral dis- 
turbance is equal to the mean velocity at the inflection point. 

5.2 Analytic solutions 

The chief analytic feature of the inviscid equation is the singularity 
at li~ c , a singularity which is not present in the Orr~Sommerfeld equation. 
This singularity is called the critical point, and its location is denoted 
by. Ajc. . Since c is complex, so is . The mean velocity V. is real, of 
course, but it may be analytically continued onto the complex plane. The 
continuation may be carried out by a power series expansion of H . 
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5.2.1 Method of Frobenius 


Two approaches have been used to find analytic solutions of the 
lnviscid equation. The first approach applies the method qf Frobenius, 
and was used by Tollmien (1929). The two solutions are . 




(5-T) 


4 f rh' PtC r>' h c i~ ■j‘' 1 

Me 


(5.8) 


where 


R Cv - | +* - 2 - C^j- ■+ (X + *gc> Z + * • * 

( 5 . 9 ) 

’ 1 *(~. * + i 4 - ■ • 

UK / 

The first solution, (p, , is regular, but because of the logarithmic term, 
is not regular near . Also, there is an ambiguity about which 

branch of the logarithm to use for aj x ^ , +TrJ or-Tri * The reso- 
lution of this ambiguity is discussed in Section 5*4. 

5-2.2 Power series in 0 1 1 

The second method of solution of (4.21) is due to Heisenberg (1924). 
In this method, the solutions are obtained as convergent power series in 
o( . These series were reintroduced into the stability problem by Lin 
(1945 ), and form the basis of his computational procedure. They have the 
further advantage of permitting simple approximate formulas to be de- 
veloped for such things as the critical Reynolds number. Their drawback 
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is that they are only convenient to use for small values of t and 
they do not show the analytic nature of the solutions as clearly as do 
Tollmen' s solutions. The two solutions are 


f.cg * (V-c'l Fi *• «! A.tcl + «! leg + • •] 

t! *- u * J 


( 5 . 10 ) 


= 0<- c ) ^i(^) + tti Ijtg) + A ijrCjV- 1 - - ' ] (5.ii) 


where 



(5.12) 




(5-13) 


= J (^77 h l 4 

The path of integration to be followed in the evaluation of these integrals 
is under the singularity, 

5»3 Numerical integration of inviscid equation 

Neither of the analytic methods described above are really adequate 
for producing numerical solutions of the inviscid equation. Only direct 
numerical integration of (4.21) can produce solutions accurately and quickly. 


There are two basic methods available. In the first, which was developed 
by Conte and Miles (l959)j> the integration is restricted to the real axis 
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and the Tallmien solutions are used to carry the solution across the 
critical point. ' In the second method, which was developed by Zaat 
{1958), the solution is produced entirely hy numerical Integration. 

This is made possible by use of an indented contour in the complex 
plane, and the calculation is performed on this contour just as easily 
as if the path of integration were the real axis. This same method will 
be applied in Part B to the integration of the compressible inviscid 
equation. 

5*4 Use of inviscid solutions in asymptotic theory 

Since Ul must be zero for an inviscid neutral disturbance, the 
second inviscid solution (5*8) is regular In this case. However, for 
a viscous neutral solution, which is constructed In part from the two 
inviscid solutions cp ( and ^ , l ( c is not zero and the logarithmic term 
of (5*8) must be dealt with. The question to answer is, which branch 
of the logarithm does one use for ^ - 1 Tr or 4 -iTT ? Tollmien 

(1929) provided the answer by obtaining a viscous correction to 
in the region around the critical point, and requiring that this viscous 
correction match fjV away from the critical point. This comparison gives 
-iff 'as the proper branch of the logarithm, which means the path of 
integration for the inviscid solution must pass under the critical point. 
The viscous correction which provides this result is derived in Section 
6 . 6 . 
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5-5 Amplified- and damped inviscid solutions as complex conjugates 
Now we continue with the inviscid theory itself. Sipce U! * o for 
a boundary layer, it follows that when > o , the critical point lies 
above the real y axis o') f and when c ; <- o , it is below the 

real axis ( ^ . < o) • When c. s o' , the singularity is on the real axis 
but thanks ’to ]l!! - 0, the logarithmic term drops out of ( 5 . 8 ) and there 
is no problem. For amplified and damped solutions, the singularity is 
off of 'the real axis, and it would seem that there is also no problem 
in these two cases. Indeed it can be seen by manipulating the inviscid 
equation, (4.2l), that if is a solution for + in; , then 

is a solution for Cy - \ C\ . Thus amplified and damped solutions 
are complex conjugates, and the existence of one implies the existence of 
the other. From this point of view the criterion for instability is that 
C is complex; the only stability is .neutral stability with c real. 
Equation (5*3) applies for c ; <0 as well as c ; > o , so unless there is 
a point of inflection, neither amplified nor damped solutions can exist. 
Since the Blasius boundary layer has no inflection point, this argument 
says that no inviscid disturbances are possible: amplified, neutral or 

damped. But viscous disturbances exist; what happens in the limit as 
R-r 0* ? ' 

5.6 Amplified and damped inviscid solutions as R -r lo limit of 
viscous solutions 

, The clarification of this point is due to Lin (l9^5)j who showed 
that if the inviscid solutions are regarded as the infinite Reynolds 
number limit of viscous solutions, a consistent inviscid theory can be 
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constructed in which damped disturbances exist and are not complex 
conjugates of amplified disturbances. To achieve this result, something 
must he surrendered, namely the existence of the damped solutions every- 
where on, the real y axis. Y/hat this amounts to is that even though 
viscosity is zero in the limit, the region over which it ' ^cts also 
goes to zero, so an integrated effect remains and the inviscid solution 
must break down somewhere. Lin's arguments were physical and heuristic, 
hut a rigorous justification was given by Wasow (1948). 

It is also possible to arrive at Lin's results from a strictly 
numerical point of view, which is the procedure to be followed here. 

To apply the method of direct numerical integration, (4.21) is replaced 
by two first-order equations for £5 (- and tU(“ (ft/ ) Tr ) . These 
equations are 

i ‘ = ~ % + — i— £* (5.110 

U-C II- c 

n * * X (5 * 15 ) 

1 " 1 ) ^3 

The analytic solutions in the free stream, 

ii = i e”' 1 '! =• -G-^e^ 3 (5.16) 

«s J 

are used to start the integration which then proceeds to the wall for 
selected values of & and c^+ic; . To obtain an inviscid solution, 
the boundary condition £ 3 C<A= O is satisfied by a linear perturbation 
of any two of & + 1 C\ . 
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The main question to decide is how to indent the contour of inte- 
gration,, The two possibilities are shown in Fig. 5*1* 



o>) 



Fig. 5*1* Alternate contours of integration for inviscid equation' 

i) For an inviscid neutral solution., cT; r o j £3 is pure imaginary 
and £4. is real. It makes no difference if the contour is indented below, 
as in (a), or above as in (b). 

ii) When the inviscid solutions are used in the computation of a 
viscous neutral solution, Cj s 0 , but £.5 Col r 0 since the boundary 
condition is to be satisfied by the inviscid solution plus a viscous 
solution. The solutions are complex and path (a) under the singularity 
must be used as shown by Heisenberg (1924) and Tollmien (1929). 

iii) If d-^|r 0 } the integration can be restricted to the real axis. 
However, then the Rayleigh theorem applies, and unless H u ~o there are 
no amplified solutions (or their complex conjugates, the damped solutions). 
But if contour (a) is used for damped solutions, and contour (b) for 
amplified solutions, both of these solutions exist even with U £ 0 . 

Some eigenvalues computed for the Blasius velocity profile are shown in 
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Table 5.1 


Contour 


Cr- 


• (a) 

0.226 

0.100 

- 0.00057 . 

00 

0.226 

0.100 

+ 0.00057 

(a) 

0*529 

0.200 

- 0.00717 

ft) 

0*529 

0.200 

+ 0.00717 

Table 5*1* 

Inviscid eigenvalues for Blasius 
with indented contours. 

velocity profile 


As can be verified from (5*1*0 and (5.15)* the solutions with c r - i c t 
and contour (a) are related to the solutions with c r +ic, and contour 
(b) by 





. Jb\ 

K * 

1 hi . 

* hr 

- i 

fis\l 

&V\ 

Jl] 


+ 


B - •? 

4 
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Which option does one pick? Since the neutral -^stability curve for the 
Blasius profile is as shown in Fig. 1.1a, the viscous solutions at all 
0< become damped in the limit R -? cq . With the inviscid solutions 
required to be the limiting viscous solutions, it is evident that contour 
(a) is to be used, just as in the asymptotic theory and in agreement with 
Lin. There are no amplified inviscid solutions. For a profile with 

:0 at , both amplified and damped solutions exist for each contour, 
unlike the Blasius case. The neutral solution is o! s / and can be obtained 
with either contour. With contour (a) the amplified solutions have 



the damped solutions. . Contour (b) gives the opposite results. 

Comparison with the viscous neutral- stability curve, which will be of 
the type shown in Fig. 1.1b, shows contour (b) must be rejected. 

The damped solutions with contour (a) do not exist everywhere on 
the real axis.' According to the theory, there is a length of real axis 
where viscosity will always have an effect and the inviscid solution does 
not apply. The numerical integration produces a perfectly reasonable 
looking solution everywhere except at , but for sufficiently large 
- c ; the numerical solution ceases to be an analytic function of • c . ' 
This matter may need further investigation as physically it is difficult 
to imagine how in the limit R ** viscosity could continue to have an 
effect over a length of the real axis and not just at a point (or several 
points). 

Figure 5*2 gives - <xc ; , the time rate of damping, as a function of 
c* s for a Blasius velocity profile. r i 



Fig. 5,2. Inviscid damping rate as function of wave 
number for Blasius velocity profile. 
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6. Asymptotic Viscous Theory 

6*1 General ideas 

We now turn our attention to finding four linearly independent 

■ 

solutions of the Orr- Sommer f eld equation at large "but finite Reynolds 
numbers, which means that an asymptotic theory in ccR is vindicated. If 
a simple expansion is made in inverse powers of t j +**), 

the first term will satisfy the inviscid equation, and thus this approach 
yields only two of the necessary four solutions. However, the inviscid 
solution designated by '$ («■ 2$} , which is the solution that behaves , 

as in the. free stream and is a combination cf and , or of 

A> " A _ 

and (J> v , goes a long way towards solving the whole problem. It 
satisfies the boundary condition at infinity; it is an adequate solution 
over the large part of the boundary layer where viscosity is unimportant; 
and the device of following an indented contour under the singularity at 
l L^C. properly accounts for the effect that the action of viscosity in 
the region of the critical point has on the solution between the critical 
point and the wall. What remains is to determine the influence of vis- 
cosity in the region immediately adjacent to the wall. A true inviscid 
solution can only satisfy one boundary condition at the wall, and indeed, 
as we have seen, only damped inviscid solutions are possible for the 
Blasius boundary layer. To satisfy both boundary conditions, we must use 
a combination of the inviscid solution which satisfies the boundary con- 
dition at infinity and a viscous correction. With the viscous correction 
designated by , and the subscript w used for conditions at the wall, 
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the boundary conditions at the wall are 


T A (^yul "** 0 
+ A <f/ w = o 


( 6 . 1 ) 


where A is a complex constant. Because these two equations are homo- 
geneous, solutions can only exist when the determinant is zero. Hence 

Mr = + (6.2) 

*-~ . §; & 

The left-hand, side is a function of o( & and c. , and the right-hand 
side is a function of dj , c j and R& . Thus (6.2) is the secular 
equation, or equation for the determination of the eigenvalues. Once 
the eigenvalues are known, the constant A is found from either one of 
the equations of (6.l), and the eigenfunctions can be computed over that 
portion of the boundary layer for which (£ and are adequate approxi- 
mations, Most of the rest of this section is devoted to finding a form 
of which will provide reliable numerical results. The viscous cor- 
rection is obtained from the simplified Orr-Sommerfeld equation 
(4.24), which was indicated in Section 4,3.3 tc be the appropriate 'iscous 
equation for large o£j, R b . This equation is, when rewritten slightly, 

cf ly - i (Ic-c) f' ~ 0> ( 6 -3) 
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The physical reasoning behind ( 6 * 3 ) is that viscosity is important 
for the disturbances only where large disturbance velocity gradients 
exist. A first-order viscous theory need only retain the largest viscous 
term (right-hand side) and largest inviscid term (left-hand side) of ' the 
complete Orr-Sommerfeld equation (^-.^ ). 

6*2 Length scales o f the viscous regions 

There are two regions of boundary layer where viscosity is of 
importance for the disturbances. ?he first of these regions is near 
the wall. Since the eigenvalues are computed from (6.2), which involves 
only wall values of {£ and ^ , the asymptotic theory is devoted pri- 
marily to developing solutions which are adequate near the wall. The 
second viscous region is near the critical point, where viscosity is 
brought into play by the singularity of the inviscid equation. As already 
mentioned, the effect of viscosity near the critical point has been in- 
cluded in the inviscid solution by the device of integrating along an 
indented contour under the singularity. However, the function does 
not correctly represent the solution in the vicinity of the critical point 
and if the critical point is close to the wall, the use of H in the eigen 
value’ equation will not be correct. The development of the viscous cor- 
rection to I near *j c is given in Section 6.6. 

It is useful to have available estimates of the thickness of the 
viscous regions. If the thickness of the wall viscous region is small 
compared with , the distance from the wall to the critical point, then 

( 6 . 3 ) simplifies to 
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< y ,v + « h$i> c - 0 

The solution of ( 6 . 4 ) which decreases with increasing is 

_c,.nfsa&, 

' = ^C«1 e 


( 6 .*) 


(6.5) 


This solution decreases to of its wall value in the distance 



We can use as an estimate of the scale of the vail viscous region. 

As decreases, it also oscillates with a "wavelength" of 2 TT £*✓ • 

A definition of the thickness of the viscous region which is more com- 
parable to the usual -definition of a boundary-layer thickness would 
give a result closer to 2 ~vr than to . The viscous solution ( 6. 5 )> 
which will reappear in Section 11.1 when we discuss the viscous insta- 
bility mechanism, is not adequate for the eigenvalue problem. 

The reason for the inadequacy of 6,5 is that the assumption c 

is too drastic. Two methods have been developed to solve (6.3) without 
making this assumption. The first solves (6.3) by the WKB method and is 
discussed in the following section. The second takes advantage of the 
fact that the velocity profile in many instances is almost linear between 
the wall and > and uses (6,3) in the form 

- ! C<j W tAj- Ajji y 1 ' = o' (6.7) 

This equation can be regarded either as the viscous equation in the 
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immediate vicinity of , or as an approximate form of (6,3) which 
is valid when the velocity profile is almost linear. It is ,: the latter 
view that is appropriate for the wall ^viscous region. The Blasius 
velocity profile constitutes a particularly favorable case for this 
assumption since is equal to zero at aj = o . With a nearly linear 
* solutions of (6.7) will he valid from the wall' to’ beyond the 
critical point. In particular* these solutions will be valid even when 
the critical point is too close to the wall for the WKB method to be 
accurate. It is the latter case which exists over much of the unstable 
region of practical importance. ' ; 

The solutions of (6,7) are functions of the variable $ * where 




Hence a new viscous length scale 


( 6 . 8 ) 


X Zu (6*9) 

has appeared* where 

e c = Rs ’wLi fe . . (6 ' 10 ' 

; , i ;i 

The behavior of the solutions of (6.7) is complicated and it is not easy 
to assign a meaningful numerical value to the constant in (6. 9). ’■ 'In- 
spection of the actual solutions* which are presented later (Figs. 6.1 
and 6.2), shows that in the wall viscous region ^3 is reduced to \/c 
of its wall value when ^ increases by 1 to 2 units of £ c „ For the 
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viscous correction to (j> near , the major .effect of viscosity 
extends over a distance of about two units of E c on each side of the 
critical point. Consequently., a length scale of 26* would be a rea~ 
sonable estimate of the half -thickness of the viscous region near , 
and also of the wall viscous region when (6.7) is the goyerning equation. 

To get an idea of typical magnitudes of £*/ and £ c and their re- 
lationship to 3 we take two points from the neutral- stability curve 
of the Blasius boundary layer: 


i) up'fer branch: - 3500* ^5 a 1«25> “ 0.37, ^c = 1# 91 

(6.1l) 

= 0.03 5jl€* = 0.098; «j c - O.I87 : 

4 ' ii \ 

ii) lower branch: = 3500., = 0.60j - 0.32,14 = 1*94, 

1 ' ; (6.12) 
£*, = 0.053,2^= 0,126, ^ = 0.162 

6.3 ~WKB solutions ; ; ;• 

The WKB method is a standard method of obtaining solutions of 
equations of the type (6.3), and was first applied to the stability 
problem by Heisenberg (1924). It will provide solutions which are valid 
far from . Let 

= (J ) . (6.13) 

and expand a ( nA 


in the series 



' f l 

Substitute (6.13 ) and (6.19) In (6.3) and equate equal powers of . 

The following two equations are obtained for and g 4 , 


8* 


IK.- A 


(6.15) 


~ (it" + ^ r 3 1 ^ 


The solutions of these equations are 


go = ± ^ 


Vz 




(6.16) 


+ X - c. 


If more terms in ^ are desired, it is necessary to return to the complete 
Orr-Sommerfeld equation, but to the above order of approximation (6.3) 
and the complete equation are the same. 

With (6.16) we find the two WKB solutions to be 


-34 


= (U' c ) e-*p [" C'tsRs)* 1 QOj)] 


(6.17) 


<$a M - (M-cS e<p[.t(«sKs) 


(6.13) 


where 


* ^ '/z 

• J [i(W-el] ^ 

V 


(6.19) 


The point , where It ~ C , is a branch point for these solutions. 
Just as for the inviscia solutions it is necessary to select the proper 
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branch by specifying the path of integration. The correct path, which 
is determined by comparison with the solution of (6. 7) to be obtained 
in Section 6.4, is the same as for t’he inviscid solution, i.e., under 
the singularity. Consequently, when c is real. 


U- C = i K~c \ f Q = | -f^v- ^ ^ 

c. * Ui-c | CJ = -Z rr ^ 


( 6 . 20 ) 


Of the two WKB solutions, ^3 decreases exponentially (and oscillates) 
with increasing ^ . Since ^ increases exponentially, it is not 
suitable for’ the boundary-layer stability problem and need not be con- 
sidered further. ^3 is a maximum at and with increasing A| it 

oscillates with rapidly decreasing magnitude. r . 1. 

With (6.15), the ratio /^(^V ^ ^ oun ^ t o be 

*/x 


, -[U s R s (U- t ')] 1 - £ 


(6.21) 


and, at ^ 5 o , 


- 


Vz - *> 


~ (tt fc Rs ^ e 


f Zt, 


( 6 . 22 ) 




The right-hand side becomes a function of a single variable if (6.22) is 
multiplied by c/llZ . 


-- sL 
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(6.23) 



The single variable' is 


I * («jR, zd* — ■ (6 - 2i| ' ) 

It iv 

It can be noted that for a velocity profile which is linear from the 
vail to the critical layer, %'w = and cjK^ - . In this case, 

i is equal to the variable -X* which is defined by (6.8). 

We can now write the boundary condition (6.2) as 

EKO - F-Ci) (6-£5) 

(6. 26) 


(6.27) 

The left-hand side of (6.25), is determined entirely by the 

inviscid solution $ . The right-hand side, F(&^ , is determined entirely 
by the viscous solution ^3 . Further, it is a universal function of i" 
and does not depend on the particular bound ary -layer profile. Equation 
(6,25) is the equation for the determination of the eigenvalues ^ , c } 

, or as it is often called, the secular equation. It is easily solvable 
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where 


t \ » hh, 


TZ f . \ ^ w 

E(d c-\ -- — — 
^ c & 


and 


F»1 



by a graphical procedure devised by Tollmien ( s 29) in which the 
functions FGl^ and are both plotted in the complex plane. 

The intersection point of the two curves determine ot if c and £ and 
hence from ( 6 . 24 ), the definition of £ , An example of this 
procedure, using the Tietjens function instead of FCjA , is given in 
Section 7 - 1 * . ' 

The only thing wrong with the method just outlined is that it gives 
poor numerical results except for £ greater than about 7* It could have 
been expected to fail on the lower branch because the wall viscous region 
occupies more of the distance between the wall and yc there, and is 
not valid near , but it is somewhat of a disappointment that it is not 
really adequate on the upper branch either until very high Reynolds numbers 
are reached. 

6.4 Hankel function solutions 

The second method of obtaining viscous solutions replaces ( 6 . 3 ) by 
(6*7) and will be valid when is a much smaller multiple of £ w than 
is permissible with the WKB method. The independent variable is $ , 
given by (6.8) and repeated here for convenience, 

5 - O3- Ri Vci^- j 6 \ (6,28) 

With the solutions designated by "X it) > (6.7) becomes 

, IV « 

a** “ 1 6 x,* - o * (6.29) 
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Pour solutions of this equation are 



(6.30) 



(6.31) 


*,« /«/**< ft Mf*] 

cCa 

-pi Vvj 


( 6 . 32 ) 


(6.33) 


where W/ } and H'ij are modified Hankel functions of order l/3 of the 
first and second' kind,, respectively. These solutions, which were first 
used in this problem by Heisenberg (1924), are usually expressed in most 
recent work (e.g., Reid (1965)) in terms of Airy functions in order to 
make use of the highly developed asymptotic theory of second -order dif- 
ferential equations. 

A tabulation of the function % and its first two derivatives is , 
given by Stuart (I963) from the paper of Holstein (1950). The real and 
imaginary parts of ^ , and also [ "X 3 1 , are plotted in Pig. 6.1 as 
functions of 5 ; the real and imaginary parts of "X5 » and | Xj i , are 
given in Fig. 6.2. The critical point is always at $=0 , and the wall 

is at some particular value of -$ . It is seen that j "X3 j , unlike 

| X3 J , does not decrease monotonically with increasing 5 . The 
asymptotic expansion of , given below by (6.3 ^)j is not valid until 













The first two solutions , (6.30) and (6.3l)> are the leading terms 
of the inviscid solutions ( 5 * 7 ) and ( 5 - 8), The two Hankel function 
solutions, (6.32) and (6.33), can he identified with the WJCB solutions 
ana g* * The asymptotic expansions of (6.32) and ( 6 . 33 ) for large 
$ are 

' x,~ ?•*?&?&'') {6 - 3k) 

- S* «p (| S 5 ' 1 e*” ) (6.35) 

If 'in the WKB solutions, ( 6 . 17 ) and (6.l8), Q(a j\ is replaced by 

Q( } 1 = « 0^)'‘|( r3 / ( 6 - 3 °) 

then 

1 P . Vi 

3 ( 6 . 37 ) 

and, since 5 - 0*j & ^0 \ ^ s e.^ , and (- i \ * e 1 '" , 

dj*5 can he identified with Xa , and <^4 with X4 • It is this identi- 
fication that establishes the integration path for the WKB solutions to 
be under the singularity, because the asymptotic expansions, (6.3*0 end 
( 6 . 35 ), of the Hankel functions are only valid in regions traversed by 
such a path, and not in the regions traversed by a path above the singu- 
larity, This whole question of the proper path to follow and the regions 
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of validity of the asymptotic expansions is highly complicated# For 
more detailed discussions,, it is necessary to consult Lin's. hook, and, 
for the complete rigorous theory, the paper by Wasow (1948)- A good 
discussion is given by Reid (l 965 ). • 

6.5 Tollmien'g improved viscous solutions . / 

The VJKB solutions are valid for ( ^ I >? > and the Hankel 

function solutions are strictly valid for = OC^A , since in 

general is not linear from - 0 to nj* . Tollmien (19^7) gave 

a pair of solutions for real c -which are valid everywhere. Similar ‘ 
solutions were obtained. for complex c by Wasow (1953)* The form of 
these solutions is suggested by a simple argument of Lin (1955)* In 
the comparison of the WKB solutions with the asymptotic expansions of 
the Kankel function solutions, the arguments of the exponentials were 
shown to be identical when iL-c in was expanded in a series about 

*jc and only the first term retained. If the two arguments are iden- 
tified without the expansion. 


? 5 /i £■"* 


| 5 “ e*'" = 


'ii. 


original page is 
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(6.38) 


then 


s * K w 


>/} 


*/j 


(lt-c\ i 


l \ 


1 


(6.39) 


With the variable 5 of (6.32) and (6,33) reinterpreted to he (6.39) 
instead o: (6.28), and with a suitable factor placed outside of the 
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integral signs,, the viscous solutions are 



These are Tollmien's solutions. Their asymptotic expansions for large 
$ are equal to the WKB solutions; for , they reduce to (6.32) 

and (6.33) (except for an unimportant constant factor). Tollmien proved 
that these solutions differ from two solutions of the Orr-Sommerfeld 
equation hy ah error of order £« for all ^ * 

6.' 6 Improved inviscid solutions 

The viscous solution discussed in Section 6.4, and , the 
solution of the inviscid equation, are adequate to solve the eigenvalue 
problem. However, for the calculation of the eigenfunctions it is nec- 
essary to get rid of the singularity in the inviscid solution at . 

This singularity does not exist in the Orr-Sommerfeld equation, and it 
can he removed from the inviscid solution by adding a viscous correction, . 

The first two of the four solutions of the simplified viscous equation 
(6.7)* which are giver by (6.30) and (6.31), are the first terms of the 
analytic series solution for and <Sf L . They are actually solutions 
of the equation <Sj '** O , which may be regarded as the inviscid equation 
with the Cp£(J and li / Ql- <■) terms neglected. Now the improved inviscid 
solution must come from a balance of the main viscous term \ (j /oij $5 , 
the convection term (It- ^ 1 } and the term ^ , which in the absence 
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of viscosity is responsible for the singularity. Hence, the equation 
for the viscous correction is 

- j ei 5 t i R 4 It ^ " o . (6.42) 

Since the scale of the viscous region in the vicinity of is 
proportional’ to £<. - (<* s f? s l/cV* s the variable £ given by ( 6 , 28 ) is 
the proper one to use* The solutions sought have the form 

■%(.?)* *?'($) + X‘"(S'i •' (6.43) 


Xth- rf’csi - 


(6.44) 


When these are substituted into (6.42) and equal powers of £ c are 

> -,£«\ 

equated, it is found that \i and satisfy the simplified viscous 

equation ( 6 . 29 ), and are given by ( 6 . 30 ) and ( 6 . 3 l). and Xi 

satisfy the equation 


tri 'V M “ 

X (i --.SV 




W 41 

■l* 


i *•* /' 

1 w *“ l 


( 6 . 45 ) 


..(rt u 

or, since A,.t = 0 > 


(il v h[ n 

- iSV 


ul 


u\ 


For Xt the proper solution is 


*V 6 ' _ l Ttr 0 

*' l ul & 


(6.46) 


(6.4?) 
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(6.46) 


since, with ( 6 . 47 ), 


x = *r + «•* x' = s * £c $* 

2Ui 


and agrees with the first two terms of the regular inviscid solution 
as given by ( 5 . 7 ). 

cn 

The second solution of (6.46), X 2 . 1 cannot be written down in a 
simple form; It is determined from (6.46) by the requirement that its 
asymptotic form for large £ behave as (* + Uj £*) , the first 

singular term of (5*8). The solution *)(J has been calculated and tabu- 
lated by Holstein (1950). The tabulation may also be found in Stuart 

<l\ Cl) 

( 1983 )- The first derivatives of "Xi*. and X^t , which are proportional 

to the real and imaginary parts of f , the amplitude function of the 

longitudinal disturbance velocity u' , are plotted in Fig. 6 . 3 . Also 

shown is the real part of the singular term of , which is also the 

asymptotic form of . It wild be recalled from Section 5-4 that the 

proper branch of the logarithm to use in for ^ <^ c could not ce 

established within the framework of the inviscid theory. The imaginary 
/n / 

part of Ax > as seen from Fig. 0 . 3 , clearly establishes the correct 
branch to be - tvi . 

According to Tollmen (1947), the _ improved second inviscid solution 
can be immediately written down as 




( 6 . 49 ) 


- 53 - 














where R and fj/ are still given by ( 5 - 9 )* However, according to 
Held (1965), Xi can only correct the first singular tepm of 

and higher -order corrections would have to be obtained tQ remove the . 

f 

higher-order singularities. Keid gives for , •' . 


t £0 )(a. ($) 


( 0 .^ 0 ) 


is close enough to the wall, so that Xi' OH differs from its 
asymptotic form . (^t /^4) J + at r o j then the viscous cor- 

rection will affect the eigenvalues, ho calculations have yet been made 

with either (6.^9) or (6.50) used in place of (5*8) for GJ Z . 

> J 

6.7 Tietjens function and the secular equation 

The ratio j ^3 > when in the form (6.32) is used for , 


in which the most important singularity is 


til 


\ 1. 



( 6 . 51 ) 


This ratio is to be evaluated at ^ -0 to provide the right-hand side 
of (6.2). In order to have a positive variable, 5 *, replaced by 
• * 

t V3 

*= - 5 * - y (6.52) 
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The ratio at *^-0 
to he. 


FCfl - 


is expressed "by Tietjens function, 'which is defined 


/Vs/V‘H;;[jwf]4s • . 

^ 32 ■ (6.53) 


cO 

Then the function which in the WKB method could he expressed as a 
universal function of i , is 


c l 1^ 

H' w y F(i) 


(6.5*0 


Define 

i * \w 



(6.55) 


If the velocity profile is linear from = o to aj. & } the quantity A 

will he zero. In any case it can he expected to he a small quantity 
and could well he ignored. The secular equation, (6.2), now takes the 
form 

kL au. = - iik 2is = dt\) F(i\ (6.56) 

c . C ‘j’j'v/ 

or ‘ "r- 

■ — EK ,0 = F(*\ (6.57) 

l + X 
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This equation is the counterpart of ( 6 . 25 ) in the WKB method. The 
right-hand side is a universal function of £ and is independent of 
the velocity profile,* and the left-hand side is a function 'of and 
c, and depends on the velocity profile. 1 

When the Tollmien viscous solution ( 6 . 4 o) is used for ft a secular 
equation of the form (6.57) is also obtained. However, the variable 
* and the function X have different interpretations. They are 

Jf 5 


% - 


% 

Ity-cl i 


1 


( 6 . 58 ) 


I + X — 


. UJ 


1 Vo 




< 6 . 59 ) 


The Tietjens function, first introduced by Tietjens (1925) in the 
form , is perhaps the most computed function in fluid mechanics. 

The paper of Holstein ( 1950 ) includes tabulations of the results of his 
own computation plus the results of five other authors. The most recent 
and extensive computation is by Miles (i960). The Tietjens function is 
equivalent to the function F(&) of the WKB method, but has an entirely 
different behavior except for large £ . The same graphical procedure 

as described in Section 6.3 can be used to solve ( 6 . 57 ). The use of the 
Tollmien solutions requires the computation of the integral appearing in 
(6.5$) and (6.59) as a supplementary calculation. 
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7- Numerical Examples of Asymptotic Method 


7-1 Calculation of neutral-stability points 

After all of the theoretical development , we can take a "brief respite 
"by working out some numerical examples. The first example will consist of 
the calculation, by Tollmien’s graphical procedure, of the upper and lower 
branch points of neutral stability for the Blasius velocity profile and ■ 

Ct- = 0. 360, The calculation proceeds as follows : 

a) From a tabulation of the Tietjens function, such as given by Miles 

(i960) or Mack (i960), is plotted against F,.(*0 to give the’ 

curve labelled Tietjens function in Fig. 7*1* 

b) From a computer program that integrates the inviscid equation, the 

function £ (et 5l c.‘)/0 + >0 of -the left-hand side of the secular equation 
(6,96.) is computed for c r = 0.3 co , * 0 and several values of <* s , 

is plotted against Ev- on the. same graph as FCs- - ) . The result is 

t 

the curve labelled inviscid in Fig. 7*1* 

c) The intersection point of the two curves establishes ot A and Z for 
the two neutral-stability points. 

The numerical results are (with $ \fR^/x= Co): 

= 0.1823, u\ = 1.922 for Cr = 0,360 

Lower Branch Upper Branch 

^ - 0.7290, %, « 2.642 = 1.204, = 3.965 

y. 

Since 2 * («{j R s l£ 1 4 ') 

«f s R» = 1584 
Rj, = 2173 
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«4Rs = 5350 
R* = 4444 








These results may "be compared with those obtained by Kaplan (1964) from 
the numerical integration of the Orr-Sommerfeld equation. 


« 0,765, R. = 2360 


*£* = 1 . 20 , = 4360 


The close agreement on the upper branch is probably a lucky accident. 

The differences on the lower branch are about what one would expect. 

If the calculations are performed with the Tollmien viscous solution 
of Seetio.* 6.5, the values of ^ remain unchanged because the left-hand 
side of the secular equation is unchanged except for the unimportant 
function X(c) ( \ changes from 1.009 to 1.0055)* Only the definition 
of 5 , and hence , changes. According to (6.57 )j 


l h 


a = K Rjl 


For 0.360, 


I* 

■> * 


Va 


l h 


i l*-‘i 4 


5 “ - Q. L1 12. fM 


V 3 


and 


- 1555 
= 2133 


^ = 5245 

Ri = 4356 


The already good agreement on the upper branch is made even better; the 
discrepancy on the lower branch becomes a bit greater. The Tollmien 
solution is an unnecessary refinement for the Blasius velocity profile, 
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since- lt ( ^ is so nearly linear. However,, it should offer a distinct 
improvement over the simpler solution for those ..profiles where U(t^) 


is far from linear. 


7*2 Thickness of viscous regions 
We can compute the viscous length scales 
£* arid £ c from their definitions (6.6) and (6.10). 


tw = 0.0592 » 0.0322 

l£ c - (2)(0. 0690) = 0.138 « (2)(0.0440) = 0.0880 

Since the functions X 3 and X t are available to us, we may examine the 
extent of the viscous regions directly from Figs. 6.1 and (6.2). The wall 
values are . 

ji/j =.7-3P • |X/j=15.^o 

t ' 

Xs is reduced to 1 ja of its wall value at 

$ = -‘ 1,0 S = - , 2.1 

Hence, in terms of ^ the thicknesses of the viscous regions are 

= (2.64 - 1.00) & = (3.96 - 2,10) £. 

= 0.113 ' = 0.082 

and 

= 0-62 flv/J, =0.45 
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These values for <^ v are not too far from the value of Zz c computed above. . 
It is apparent that £ w gives a considerable underestimate of the thickness. 

However, it should be noted that the ratio of at the two neutral points 

* 

is proportional to neither the ratio of £ w or . 

- By reference to Fig. 6 . 3 for the correction function Xj* * > we* are 
able to check directly whether the critical -point viscous region extends 
to the wall. For the upper-branch point {z - 3* 96), it is evident that 
this region does not extend to the wall, although it does overlap the wail 
viscous region. However, for the lower-branch point ( - 2.64-), *X iv% is 
still 2 jfp "below its asymptotic value at ^ s 0 • It is not known how much 
of an effect this small difference would have on the eigenvalues, but we 
can -expect any effect to be well within the order of accuracy of the method. 
7*3 Calculation of an eigenfunction 

The final step in the computation is to obtain the eigenfunctions, and 

1 

we ohall partially indicate how this is done for the longitudinal velocity 
fluctuation at .he upper -branch point calculated above. The inviscid ampli- 
tude function F- < is plotted in Fig. 7*2. It is valid for y greater 

than + d£; (3t* = 0.132). For ^ c ~3£« t t 4 H« it must be 

corrected by Xj J,/ , and for ^ ^ ^ ft must also be corrected by X 3 . 
Consequently, for ^ ^ these two corrections must be applied simultane- 
ously. 

The first step is to compute the complex constant A (actually iA here) 

' . ji’"' 

from the second equation of ( 6 . l). * - - 
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We obtain Fv and F; from the output of the computer program for 
C* = 0.360, ~ 0, ^ « 1,204 The ’X^ and X 3 . are taken from 

i 

Stuart's ( 1963 ) table, with the second derivatives used to calculate 
the values listed here. 

F =-5*217 • Fr^ =-1.915 

y'co = 6.09 y ' <m = - ik.ik 

With these values, we find 

A r = 0.0198, A*, = 0.360 

The corrected amplitude function is found from 

■fc^ = Fty * AXj'tjJ 

and is shown in Fig. 7*2. 

In Fig, 7*2 the viscous correction in the critical -point region is 
applied only to the imaginary part of F * This is done quite simply, 
hy multiplying o4* by a constant C such that C Xi| d. 131 ^ f 
and then continuing F; heyond ^ -0-132. with C * To correct 

Fv it is first necessary to separate 1 + t&j. $ from F^ and this has 
not been done, A complete description of the correction procedure may 
he found in Schlichting (1935)* where the inviscid solution through the 
entire boundary layer is obtained from the series expansions in powers of 

V’-K" 
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8. Numerical Integration of Orr-Sommerfeld Equation 

The asymptotic method described in the previous section provides 

eigenvalues of neutral disturbances with reasonable accuracy. However, 

* 

for the best results it is necessary to have available a computer program 

for calculating exact numerical inviscid solutions, and, in the absence 

of computations of the Tietjens function for complex S , simplified and • 

damped solutions must be obtained by a different- method which involves 

a perturbation of the neutral solutions (Schlichting ( 1933b), She.i (195*0)' 

Because of these limitations, attention has turned in the last few years 

to obtaining all numerical results from exact numerical solutions of the 

Orr-Sommerfeld equation itself. The task of developing suitable numerical 

methods was not an easy or', because in the Orr-Sommerfeld equation the 

highest-order derivative ^ _s multiplied by a small quantity, * 

This circumstance, which is the source of the non-regular nature of the 
/ 

inviscid and \ IKS solutions, is also a source of great difficulties in the 
numerical integration. 

The WKB solutions clearly demonstrate the reason for this difficulty. 
It can be seen from the solutions and cf* , given by (6.17) and (6.18), 
that regardless of whether the integration proceeds in the + y or - y 
direction, there will always be a solution present which grows roughly as 
& ] , Any eigensolution of the Orr-Sommerfeld equation is 

of order one throughout the boundary layer, and it is the task of the 
numerical integration to produce this solution in the presence of’ solutions 
<$3 and <jV . Since every numerical integration scheme has some inherent 
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errors, and the growth of these errors will he determined by the most 
rapidly growing solution, it can easily happen that after the inte- 
gration proceeds a certain distance, .the error will overwhelm the 
solution. 


The available numerical integration methods fall into three main 
types. The earliest in time (for channel flow) was by Thomas (1953 ) f 
where the forward integration method was abandoned at the outset. Instead, 
a finite -difference approach was used where the Orr-Sommerfeld equation 
was simultaneously solved at a large number of y stations, and the eigen- 
values obtained by the inversion of a large matrix. This method was 
later successfully applied by Kurtz and Crandall (1962) to the Blasius 
boundary layer. The second approach is that of Brown (1954, 1959)* He 
used a straightforward forward integration method and errors were con- 
trolled by restriction to low values of ^ i? s , or by an increase in 
the number of significant figures (use of double -precision arithmetic). 

The third approach was initiated by Kaplan (1964). In this method, the 


integration starts in the free stream and proceeds to = o . Two lin- 
early independent solutions are produced. The first solution, which for 
y h is equal to the viscous solution (4.19), integrated first. It 
becomes large and offers no problem. The second solution, which for 
? I is equal to the inviscid solution; (4.l8), grows much less than 
the first solution and is in danger of being "contaminated" by it. To 


avoid this, a quantity proportional to the first solution is subtracted 


from the second solution after a few integration steps, and a condition 
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is applied which forces the second solution to remain small. This 
process was called "purification" by Kaplan, and a version of the 
method is described in detail in Betchov and Criminale (1967). A 
rigorous version of Kaplan's method was provided by Bellman and Kalaba 
(1965) and applied to the boundary layer by Wazz&n, Okamura, and Smith 
(1966). In this method, the standard Gram-Schmidt orthonormalization 
technique is ‘used to assure that the second solution is linearly in- 
dependent of the rapidly growing solution and of magnitude unity. 

Another method along these lines that is worthy of mention is the 
quasilinearization method of Raabill and Van Driest (1966 ). 

The method of Kurtz and Crandall and the forward integration method 
with single -precis ion arithmetic (8 digits) are limited to growth factors 

g 

of the rapidly growing solution of about 10 . Double precision in the 

l8 

forward integration method extends this growth factor to about 10 

In contrast, according to Betchov and Criminale {1967), the growth factor 

IfO 

in .Kaplan's method can .he 10 , and in Wazzan, Okamura, and Smith's method 

80 

it can be 10 . Both of these methods use only single precision arithmetic 

The latter method can be considered to solve the problem of integrating 
the Orr^Sommerfeld equation under almost any conceivable circumstances. 

For example, the latter authors have computed the eigenvalues of neutral 
and amplified disturbance? for the Blasius boundary layer up to R* - 
where ft*, is the x Reynolds number. Instability first appears at f'q; Z. f* io*j 
and transition, even for a free stream of very low turbulence level, occurs 
by R* a 3 * io c . Therefore, this method more than covers the range of 
Reynolds numbers of interest for the transition problem, and could be used 
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to study such things as the approach of finite -Reynolds number damped 

solutions to the.inviscid limit with increasing Reynolds number. 

No details of any of these methods will be given here.! In Section 

lb, 2, a brief description is presented of a forward integration method 

r 

which provides automatic calculation of the eigenvalues and eigen- 
functions of the compressible stability equations. 



9 « Numerical Results 

9.1 Reyiev of existing calculations 

A- co.nsid.er able body of numerical results has been ‘accumulated over 
the years from both the asymptotic method and from direct numerical inte- 
gration of the Orr-Sommerfeld equation. These results, while extensive, 
did not begin to be comprehensive until the computer came into use. The 
earlier investigators usually contented themselves with the computation 
of the neutral- stability curve. A notable exception was the calculation 
of Pretsch (19*0.), where amplified as well as neutral solutions were ob- 
tained for a family of velocity profiles. After the development of 'approxi- 
mate formulas by Lin ( 19 * 1-5 ), it became the somewhat unfortunate custom to 
compute only the critical Reynolds number, i.e., the lowest Reynolds number 
for instability. Such 4 calculations, as has been emphasized by Kaplan (196*+), 
are of limited usefulness. In order to have an accurate’ picture of the 
comparative instability of two different boundary-layer profiles,, it is 
necessary to calculate the neutral -stability curve and the amplification 
rates. In the following summary, only those calculations are mentioned 
in which at least a neutral-stability curve has been obtained. 

The neutral -stability curve of the Blasius boundary layer has been 
computed, using the asymptotic method, by Tollmien (1929), Schlichting 
(l 933 a ).> Pretsch (19*0.), Lin ( 19*^5 Lees ( 19 **-?)* and- Shen ( 195*0 J and, 
using numerical integration, by Brown ( 1959 ), Kurtz and Crandall (1962), 
Kaplan (196*1), Betchov (1965), Vazzan, Okamura and Smith (1966), and 
Radbill and Van Driest (1966 ). The effect of pressure gradient was first 
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considered "by Schlichting (1940) for the Howarth profiles . 

Neutral-stability curves for the Falkner-Skan family, V.t ~ X , have 

\ 

been computed by Pretsch ( 194-1 ), Kaplan (1964), and Wazzan, Okamura 
and Smith (1966). Neutral- stability curve., for 6th-degree polynomials 
have been computed for various values of the Pohlhausen parameter by 

1 

Schlichting and Ulrich (1942) and by limman, Zaat and Burgerhout (1956). 
The boundary layer with suction was computed by Pretsch (1942), Bussmann 
and Miinz (1942), Freeman (unpublished, but see. Chiarulli and Freeman 
(1948)), and Hughes and Reid (1965). Contour lines of constant non- 
zero C] . from which the amplification rate may be obtained, have been 
computed for various flows from the asymptotic theory by Schlichting 
(1933b) and Shen (1954) for' the Blasius boundary layer, and by Pretsch 
(l94l) for the Falkner-Skan family of profiles. Such calculations are 
usually performed as a matter of course when the method of numerical 
integration is used. Eigenfunctions have been computed by Schlichting 
(1935) and Holstein (195°) from the asymptotic method, and are normally 
obtained along with the eigenvalues in the numerical integration methods. 

Some results for the Blasius boundary layer will be found in Section 
12, where the experimental confirmation of the theory is presented. The 
calculations performed for the Falkner-Skan family of velocity profiles 
are of particular interest because these profiles provide examples where 
1 

Doubt has recently been expressed about the validity of calculations 
made from the parallel-flow theory for the pas.e of a favorable pressure 
gradient by Morkovin and Donaldson (private communications). 


inviscid instability; is important, and thus the results serve as 
background for the stability of the compressible boundary layer to 

be presented later. A short discussion of these results is given in 

r . 

the next section. 

9.2 Falkner-Skan velocity profiles 

The Falkner-Skan family of boundary-layer velocity profiles cor- 
respond to free-stream velocity distributions of the form 

fc, 00 * k X* 1 (9-1) 

The boundary -layer equations reduce to ordinary differential equations 
for these profiles. These equations were solved for various values of 
•m by Hartree ( 1937 ), and more exactly and with more detailed results 
by Smith ( 195*0 • The similarity variable which replaces vj. is 



(9-2) 


where 

0 = 

h.-1-l (9-3) 

The value ^ r o corresponds to the Blasius boundary layer; with ^ > 0 
there is a favorable pressure gradient; with ^ < 0 there is an adverse 
pressure gradient. The value ^ corresponds to the separation 

profile, where U' w - 0 
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According to the theorems presented in Section 5*1.> "the necessary 
and sufficient condition for neutral and amplified inviscid solutions 
is that lL U must be zero somewhere in the boundary layer. For 0 , 
lc‘'*Q everywhere (except at for the Blasius boundary layer). 

However , for ^ < o } V will always be zero at some. ^ . Consequently, 

for 1^40 the Falkner-Skan profiles have inviscid instability and the 
curves of neutral stability will be of the type shown in Fig, 1,1b. 

These profiles are of interest because they provide an example where 
there is a steady evolution from a profile where viscous instability • 
is dominant ( ^ - 0 ) to profiles where inviscid instability becomes 

increasingly important. In Part 3,, we shall find that a similar evolution 
takes place with increasing Mach number for the compressible boundary 
layer, on an insulated flat plate. 

Table 9*1 gives some numerical results for the Falkner-Skan profiles 
' * • 1 

with adverse pressure gradient. In the table, is the similarity 

variable as defined by (9*2) evaluated at s , which is 


$ 

V * 

Cs 



Vi 


0 

0 

0 

0 

4.4 x 10" 3 

6.0 

1-7 

- 0.05 

- 

0.32 

0,30 

5.9 x 10" 3 

- 

1.8 

- 0.10 

1.3 

0.43 

0.46 

1.0 x 10~ 1 2 

6,5 

2.0 

- o.i4 

2.1 

0.49 

0. 64 

1.6 x 10 " fi 

6.7 

2.2 

- 0.1988 

3-3 

0.50 

1,31 ' 

6.3 X 10~ 2 

7.8 

3.2 


Table 9*1* Numerical results for Falkner-Skan velocity profiles. 
Adverse pressure gradient. 

1 

These results, and Fig. 9«lj were obtained from material kindly 

furnished to me by Mr. A. M. 0. Smith of Me Donnell -Douglas Corn, Long 
Beach, Calif, and computed by Wazzan, Okamura, and Smith (i960). 
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the ry where = 0 . The third column,, Cs , is the phase velocity 

of the inviscid neutral disturbance, and is equal to the mean velocity 

• * 

at r^ s as was shown in Section 5*1* The fourth column gives the 
corresponding wave number (<*$*)? referenced here to displacement thick- 
ness, as obtained from the upper branch of the neutral -stability curve 
at high hut finite Reynolds number. It could be computed exactly from 
(4.2l) for c r 's c s , C;*. o ' . ' The fifth column of Table 9*1 gives 
the maximum time rate of amplification ■ (at finite Reynolds number except . 
for ^ = - 0.1988), again referenced to displacement thickness. The final 
two columns give the boundary-layer thickness and the displacement thickness 
both in the ^ coordinate. 

It is seen from Table 9*1 that with increasing negative values of 
^ , the quantities , c s , and (<=< s » c 5 ) M(a increase. Thus, insta- 

bility increases as the inflection point moves away from the wall. The 
amplification rates of Table 9*1 for contrast strongly with the maximum 

amplification rate of the Blasius boundary layer. It is evident that a 
velocity profile with an inflection point located well away from the wall 
has much stronger instability than a profile with no inflection point. The 
most unstable profile of all is the separation profile. 

The neutral- stability curves and contours of ^ - const, are of great 
interest. ' Figures 9-1 and 9.2 give these results for ^ - - O.l^f and 
- 0,1988 (separation profile), respectively. Figure 9*1 clearly shows 
that in spite of the large inviscid instability, there is still a small 
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destabilizing effect, of viscosity. . The maximum amplification rate 
at a fixed increases about 15$ as Rf decreases from infinity* 
to about 10 . Note that the maximum amplification rate is ‘-nearly ■ ' 

‘ k ' 

constant until Rs* is reduced below 10 . Most of the C\ £ const. 

v * 

c 

contours are open for ^ ~ oaA- 9 i.e., the upper and lowqcr branches 

.extend to infinity. Only those contours with c^o.oa- fqrm closed 
loops. .In contrast, for the Blasius profile (^ = 0 ) all of the contours 
except for CjsO close at seme finite Reynolds number. 

Figure 9*2, for the separation profile, shows the final stage in 
the evolution of the - const, contours. Here none of the contours 
are closed, and they all extend to infinity. Also the upper branches 
are all either parallel to the axis or have a positive slope. As 

a result, the maximum amplification rate at R-v & is the highest ampli- 
fication rate, that' is possible for this velocity profile. The destab- 
ilizing ’effect of viscosity has completely disappeared, and viscosity 
serves only to damp the disturbances. These disturbances are formed 
by the inviscid instability mechanism associated with the inflection point. 
Neutral -stability curves of this type are common for free -shear flows, 
where the inviscid amplification rates are even larger than for the 
separation profile. We shall find in Part B that they are characteristic 
of the compressible boundary layer at a sufficiently high Mach number even 
in the absence of a pressure gradient. ■ '■ 


We use the term "destabilizing effect of viscosity" to mean that the 
maximum (with respect to at a fixed ) amplification rate increases 
with decreasing Reynolds number over some portion of the Reynolds number 
range . 1 v\> * v - 
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3.0. Kinetic -Energy Equation 

Some of the e'arliest attempts at a stability theory were carried 
out on the basis of energy methods. These methods vere aimed primarily 
at a calculation of the critical Reynolds number, but vere not very 
successful even in this limited objective, as the Reynolds .numbers 
obtained vere always extremely lov. This was a , direct result of the 
fact that the disturbances used in these analyses vere not required to 
be solutions of the differential equations. Consequently, the con- 
ditions necessary to secure stability of all disturbances had to be 
more severe than those necessary to insure the stability of only hydro- 
dynamically possible disturbances. For our purposes the energy method 
will be useful in shedding some light on the physical basis of the 
stability theory. ' Up to nov, the discussion has. been mainly mathe- 
matical, and the time has come to bring out some of the physical 
mechanisms at work. ' 

This section is devoted to a derivation of an equation for the 
kinetic energy of a disturbance. The dimensional x-momentum equation 
is, in two dimensions, ■ 


„ In . l Op 

+ ■ U + V 5 = - J~ -J- -r If 

u ** h e 3 * 


/Yu d\ \ 
Hi 


( 10 . 1 ) 


When the flow quantities are split into a steady mean component and a 
disturbance, the resulting equation is * . 


V 

U 


— - + 

u' 


*. v' ii 1 V' . % 




9 7 

. Vel - 

f 

IP 

/Vu' A' . 3'fc-' 

■f- ir 4- * — 

3X 


ox 

U* 1 y y. 




( 10 . 2 ) 
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Three mean-flow terms have been dropped by the 'boundary-lay^r approxi- 
mation, but nonlinear terms have been retained. These tenps can be 
transformed as follows with the use of the continuity equation. 

* l 

* 

u'2t' * V< = — (uV) + ; ' (10.3) 

u ^ sx\ $ 

The disturbance is assumed to be periodic in X and growing in time. 

Thus an x average taken over a wavelength at a particular time has 
a definite meaning. Since such an average of periodic terms is zero, 
only nonlinear and mean-flow terms remain when the equation is averaged. 

t • 

■?/ 9k lr ^ . Q /T"7i l <?P $k ' fin l.\ 

It + V + V I „ ; — •+• IT z - Uv.t; 

2 % ^ ^ ( 9 * 3 ^ 

The quantity . 

* „ * - ( 10 . 5 ) 

is the Reynolds stress. It is a momentum transfer caused by the dis- 
turbance, and acts in addition to the mean-flow stress ^ alt/3^. If 
( 10 . 4 ) is written in the form 


p u + 
V 3 X 


pV 3 



( 10 . 6 ) 


we recognize it as the mean flow boundary-layer equation with an additional 
stress 'term. Consequently, when nonlinear terms are included, the mean 
flow in the presence of a periodic disturbance cannot be the same as the 
mean flow with no disturbance (see Stuart (1956)). 
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The next step is to subtract the mean -flow equation from (10.2), 
and then to multiply the remaining equation by u' in order to form 
part of a hinetic energy equation. This equation is 



-t* 


m' JL (u.'v*) + 



U # JL fiTT 7 ) s + ir 


(10.7) 


When the same procedure is applied to the y momentum equation, except 
for multiplication by \r' instead of u' , and the resulting equation is 
added to (10.7)^ we have 




< t ,*• 
n* v 


u'v' 





+• u 


2 fwV ) +v'j>(uV! 
•7 3 


+ u’ 


'Hi 

n 


(t 


.< -av 1 


1 (uV) 



( 10 . 8 ) 


j_ A*' v 1 hi \ + s- v' yV' ) 


The last two terms of the first line can be converted into energy fluxes. 
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The energy equation then becomes 




4 * H 


L\ (S* ?') * ^L(^) 

3jC > V 2- / 3X \ t J 


*• JL (£*jd ) 

3 J V 2* ' 


' 3v'* 4 . u'v' 


- v IT + *’V __ 

3^ <?*£• 3 


(lO.lO) 


i (u 1 It' + \r' JZ''\ + ^(n'^V + \r‘ V^v’ ) 

3 X I 


At this point we integrate ( 10 . 10 ) over a volume which extends 
one wavelength X in the x direction and to infinity in the y direction. 
The terms with a single fluctuation quantity, or an x derivative of a 
fluctuation times a mean-flow quantity, are zero by periodicity. There 
are two terms of the type u'36/3* + v l Kh^ f Their integrals can be shown 
to be zero as follows : . 



From periodicity and the boundary condition at ^ -**■ « v' , we have 



(10.12) 


This integrand on the right-hand side is zero by the continuity equation. 
The viscous terms of (10.10) can be transformed into 

v'll' •• • . 

3^, 9X (10.13) 


VX/ +■ U / 
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vhere 


ki _ ^ 21* 

3^ ~ Tit 

is the disturbance vorticity. 
integration by parts, 

' „ v ( 2^ Wx4*. - 

3 ^ U I 4 



( 10 . 14 ) 


The integral of (10.13) gives, after 





( 10 , 15 ) 


An alternate form of this integral, which is perhaps more easily tdenti 
fiable as the viscous dissipation, is 



(10.16) 


With these reductions, the remaining terms in the energy equation are 


Now, 


x* \ 

-nJ " v ' u - -? 7 


V*' a -V 


(lO.lB) 


where the bar refers to an x average, and 'V has already been identified 
as the Reynolds stress. Hence, ( 10 , 17 ) can be written, with asterisks 
introduced again to identify dimensional quantities. 
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where 



(10.19) 


( 10 . 20 ) 


is the distur bance Kinetic energy in the volume of fluid. 


VJ 



is the energy production term, and 





( 10 . 21 ) 


( 10 . 22 ) 


is the energy dissipation term. 

If .all quantities are made dimensionless in the usual manner, with 
lii as the velocity scale, and S as the length scale, then the integrated 
Kinetic -energy equation becomes 


■ _ p _ J- V 

where 



©■ -- s-y 


( 10 . 23 ) 


( 10 . 24 ) 
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and 


T? « - iaV 


(Xo/25) 


are the dimensionless counterparts of (10.20), (10.21); (10.22) and 


' ( 10 . 18 ). 

Equation (10.23) is exact in the sense that no further approximations 
were made "beyond the boundary-layer approximation. It holds for nonlinear 
as well as linear disturbances; although; as already mentioned; In the 
latter case U must be interpreted as the mean velocity in the presence 
of a disturbance and differs from the disturbance -free mean. flow. For 
a linear disturbance; as seen from (10. 6); U is the same with a dis- 
turbance as without one. The nonlinear and pres sure -gradient terms of 
(10.10) contribute nothing to the integrated energy balance. They serve 
only to redistribute the energy within the boundary layer. 

The term F of (10.23) represents the interaction of the mean flow 
with the disturbance. If the disturbance results in a Reynolds stress 
such that "P > 0 ; then energy will pass from the mean flow into the dis- 
turbance. The dissipation term represents Kinetic energy lost by the 
irreversible viscous worK’ and is necessarily positive. Consequently; the 
conditions for damped; neutral; and amplified disturbances are; respec- 
tively* 


-P 




R, 


V 


P«- V 




(10.26) 
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The purpose here is not to construct a stability theory based on 
(10.23). The eigen solutions of the Orr-Sommerf eld equation already 
provide complete information about th*e stability characteristics. 

Rather, (IO. 23 ) can serve to point out t v e direction to go in increasing 
our understanding of the physical mechanisms which are at work to produce 
instability. It is clear from (10.23) that a decisive role is played by 
the Reynolds stress and it is to this quantity that we must turn our 
attention* 
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11. Distribution of Reynolds Stress 


11.1 Prandtl ( s analysis of vail region 

Prior to 1921 it bad been accepted doctrine that a flow which was 

■o 

stable in the absence of viscosity could only be made more stable by ' 
the presence of viscosity. Although an earlier,, and in this connection 
evidently unnoticed, paper by Taylor (1915) pointed out the possible 
destabilizing influence of viscosity, it remained for Prandtl (1921) 
to clearly demonstrate the mechanism by which a stable invite id flow 
can be made unstable by viscosity. It was this discovery which led to 
the elaboration of the stability theory in later years by Prandtl 
students and colleagues. In view of the importance of this discovery, 
it is of interest to quote a few lines from Prandtl ’ s paper (my -trans- 
lation ) . 

"Previous mathematical investigations on the origin of turbulence 

t * * 

have led to the opinion that small disturbances of a viscous, laminar 
flow between two walls are always damped... In order to learn how i tur- 
bulence actually originates, I bad built at Gottingen an open channel, . . 
and observed the flow by the Ahlborn method (sprinkled lycopodium- powder) 
...wave forms with slowly increasing amplitude were occasionally observed 
...these waves of increasing amplitude completely contradicted the dogma 
of the stability of laminar motion with respect to small disturbances, 
so' that at first I tended to believe that I had not seen this infrequent 
phenomenon completely right." 
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"We now applied out selves to the theoretical treatment, and, to 
anticipate a little, we found, contrary to the dogma, an instability 
of the small disturbances. " 

Prandtl’s argument was later repeated in a different form by Lin 
(195^ 1955 )j but we shall follow essentially the original derivation 
here. The viewpoint adopted here is the same as that used in Section 
6.2 to derive viscous solutions near the wall, i.e., it is assumed that 
h(^\z 0 „ With this assumption, the x -momentum equation, from (3*l)j 


simplifies to 

3 * 7 OX 


1 .J 

(11.1) 

1. .1 


where the term v'*U/4^ has also been dropped, 
viscous region, (ll.l) reduces further to 

w ' 

‘U D* 


Outside of the wall 

1 i 

(11.2) 


The velocity w' consists of two parts; an inviscid part u- which 
satisfies (ll. 2), and a viscous part which .satisfies the difference 
between *(ll.l) and (11.2) and together with the no -slip. boundary 

condition. 

0 u J t ^ kv 

u = v (11-3) 


The solution of (11. 3) is, for = 0 and hence .<^ real. 


w v ( ^ 


_u : coi e 


(>“ 1 1 If ^ 


Cll .4) 


-O5- 



where the boundary conditions that mV) * u/ta} + ut t °\~o and u'(^\ 

outside of the viscous layer have been applied. 

The additional longitudinal disturbance velocity u f v , which is 

* 

needed to satisfy the no-slip condition, induces an additional normal 
disturbance velocity \rl through the continuity equation, 


v v C^\ = ~ 


2Uv j 

* 


(11.5) 


which yields, upon substitution of (11.4), 


v; 


,‘(^1 = Cl-Htfj «;(»!- 


(ifaRi) 


l ft 




"] 


i C * jX - ^ 


( 11 . 6 ) 


Outside of the viscous region vj is independent of y 


, c(j Ui M ; t £) 

v v - — c * — * j 7 e 




% 


' (ll-T) 


and is zero. 

The consequences of (ll,7) for the Reynolds stress, are as follows. 

We have .noted in Section 10 that the Reynolds stress is given by = - u'v-’ 

(in dimensionless form). For an inviscid neutral disturbance ana v'. 
o 

are 90 out of phase and is zero. However, for any othei disturbance 
U* and v* are correlated and there is a Reynolds stress. Since u(, is 
zero outside of the wall viscous layer, it can .contribute nothing to t; 
outside of the layer. However, vj persists for some distance outside 
of the wall layer, and, since it is shifted in phase 135 ° with respect 
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to u{ , it will ' produce, a Reynolds stress. This Reynolds stress must 

equal the Reynolds stress set up Toy the flow in' the vicinity of the 

i . 

! critical layer, and which, in the absence of the wall viscous region, 

[ .• i 

| would extend to the wall. This stress is discussed in Section 11.2. 

The formula for the -Reynolds stress at the edge of the wall viscous 

region can be derived from (11.7)* ‘If we recall that it is necessary to 

take the real parts of and vj before multiplying and averaging, we 

find 

| ■ ‘ 

= T-^-T Kw] 1 (11-8) 

l 

where - ^tcy. if the ratio is formed, we have 



(H.9) 


f ■ ... i . 

A general expression for ^ in the wall viscous region can be obtained 

] :: : 

from (il.4) and (ll,6), and this expression would give the increase of 

- • £ i- 

0/ from zero at the wall to the value given by (ll,9) at the edge of 
the viscous region. Equation ( 11 . 9 ) establishes that *€ is positive 
and, according to (10.2), energy will be transferred from the mean flow - 
to the disturbance. Consequently, the wall viscous region, which is 
formed to satisfy the no-slip boundary condition, also has the effect 
of creating a Reynolds stress which acts to destabilize the flow. This 
mechanism must be present to some extent for all disturbances, but whether 
a. particular disturbance is actually amplified; or damped will depend on 
the distribution of the Reynolds stress through the entire boundary layer 
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1; 


and on the dissipation term. 

As a note of caution, it must he recalled that the preceding 

m 

analysis rests on the neglect of li in the vail viscous region. 
Consequently, ve can expect the results to he valid only at high 
values of Rs , vhere the vail viscous region will he thin compared 
to the boundary -layer thickness, and for values of ^ large compared 
to the thickness of the wall viscous region, 

11.2 Reynolds stress distribution from the irviscid theory 
The inviscid theory can also provide us with a formula for the 
Reynolds stress which will apply to all of the houndary layer except 
the' two viscous regions. First, we need a general expression for the 
'Reynolds stress. Since ■ ■ ; . :■ 



• «(«{* <“ A'*) 

*? <2 


/ « i ( X " ^ 

V J 6 . 


( 11 . 10 ) 


ve find 

V - 


(f 


f'- $!■ 




2,<*l t 


( 11 . 11 ) 


It is easy to verii’y that the quantity V defined hy (5-4) is equal 
to the part of (ll.ll) which is in parentheses. 


9 4 y') = f- ji < u - 12 ) 

Therefore, the Reynolds stress can he written 
2 ^ ^5, c\ £ 

= i <s Ve s ' . . (11.13) 
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We have, from ( 5 - 2 ), 
C\ li 


i$r 


{U~t\ 

and it follows from (11.10) that 

V - -£ fjl e 

Consequently,, the' derivative of the Reynolds stress is 




.it 


Ci U 


( 11 . 14 ) 


( 11 . 15 ) 


This equation was first derived by Foote and Lin (1950). 

‘For a neutral disturbance, 4 ?/^ is zero everywhere except at 
U-c where, by the same derivation that led to (5.6), t has a discon- 
tinuity given by 


(11.16) 


+■ <? ] - X C nj - a ) v* L 

U'c 


(11.17) 


We can now sketch the expected Reynolds stress distribution in the boundary 
layer for a neutral viscous disturbance. 



Fig. 11.1 
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The Reynolds stress builds up in the vail viscous region to the constant 
value given by (ll.9) which, for consistency, must be equal to the value 
given by (11.17) with action of viscosity near 

smooths out the discontinuity in V , The integral of X All/ Arj must be 
equal to 3)/Rg, the integral of the dissipation terns. 

The results obtained in Section 5*1 concerning the importance of an 
inflection point in the inviscid theory also follow directly from (11.16), 
and (11.17). Since there is no dissipation in an inviscid fluid, .the 
Reynolds stress must be zero throughout the boundary layer for a neutral 
inviscid disturbance. Equation (H.17) shows? immediately that Ui must 
be zero. For an amplified disturbance, integration of (11. 16) 'from aj » 0 
to and application of the boundary conditions on yields the 

- t\ 

result that «- must be zero somewhere in, the": boundary layer. " 


:: 3 


C LI 
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12. Experimental Confirmation of Stability Theory - Schubauer- 

Skramstad Experiment 

Jfo account of the laminar stability theory is complete without 
a mention of the classic experiment of Schubauer- Skramstad (194-7). 

Up to the time this experiment was performed there was ’ ■ experimental 
evidence to confirm Tollmien’s theory of the stability of the Blasius 
boundary layer. Indeed what evidence there was indicated that the size 
of the disturbance was the controlling factor. The general attitude 
towards the theory was one of skeptic ism, and most people probably 
agreed with Taylor (1938) in thinking that the theory had little or 
no connection with boundary-layer transition. Only in Germany does 
there seem to have been general acceptance of the theory. 

This picture was radically changed by the Schubauer -Skramstad 
experiment which provided a complete demonstration of the validity of 
the theory, and also firmly established the connection of the amplified 
oscillations with transition. This demonstration was unmistakable in 
spite of the fact that the calculations in existence at that time were 
not very accurate. The exact numerical calculations of recent years 
have only strengthened the demonstration. The experiment was performed 
during 194-0 and 194-1, but' security restrictions prevented its publication 
in the open literature until 194-7 . 

Experimental techniques will be discussed elsewhere in this course 
and We are only interested here in the results. It will suffice to say 



that the disturbances were produced artificially in the boundary layer 
by means of a vibrating ribbon, which thus, provided a means of con- 
trolling the amplitude- and frequency, of the disturbances. A hot-wire 
anemometer was used to follow the disturbance as it travelled down- • , 
stream. The amplitude , and phase of u' were the primary quantities 
measured. 

The three main results of the experiments were; (i) the neutral - 
stahility curve,- (ii) the amplification rate as a function of wave 
number at three fixed Reynolds numbers, and (iii) the distribution with 
of the amplitude of w' . The comparison of the experimental neutral- 
stability points with the neutral- stability curve computed by Kurtz 
(19&) is shown in Fig. 12.1. The agreement is very good except near 
the critical Reynolds number where there is a great deal of scatter. 

In the pomparison of the amplification rates, it has for some 
reason -been customary to plot c ( rather than o£c,* . Figure 12.2 gives 
c : vs. c^*. at =12,00 as measured and as computed by Kaplan. 

What was actually measured is the spatial amplification rate, which is 
transfopned into a temporal amplification rate by (2.26). Schubauer- 
Skramstad used the phase velocity instead of the group velocity and 
achieved the remarkable agreement shown in Fig. 12. 2. The agreement at 
the other two Reynolds numbers is not as good, and at the lowest Reynolds 
number there is a factor of two difference in the maximum amplification 
rate between theory and experiment (use of the group velocity would in 
every case worsen the agreement). - ^ 
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The third comparison is given in Fig. 12.3 where experimental 
measurement of the amplitude of u' as a function of ^ is plotted 
along with an eigenfunction calculated "by Kaplan. Again the agreement 


is remarkable. 












PART B. COMPRESSIBLE FLOW 


13. Derivation of Equations 

13 • I Introductory remarks 

The theory of the stability of a compressible boundary layer differs 
sufficiently from the incompressible boundary layer to warrant being treated 
as an entirely separate subject. Many of the general ideas are the same, 
and for this reason the incompressible theory can be regarded as an indis- 
pensable prelude, but there are at least three main differences. First, 
the asymptotic theory, which served to give an adequate account for incom- 
pressible flow, has had only a limited success for supersonic flow. Almost 
all results will have to be derived from numerical integration of the sta- 
bility equations. Second, the inviscid theory has a much greater relevance 
than it had for incompressible flow, where it served mainly as a source of 
two solutions for the asymptotic theory. Although we shall concentrate 
entirely on the zero pressure-gradient boundary layer, the inviscid theory 
even here will be of great usefulness. Third, the existence’ of regions of 
supersonic flow relative to the phase velocity of the disturbance permits 
waves of an entirely different type than those found in incompi essible flow. 

The compressible theory has been developed by fewer people than was 
the case for the incompressible theory.. Aside from an early investigation 
by Kuchemann, the theory can be considered to have started with a report 
by Lees and Lin (19^6). In this work, the inviscid theory was considered 
in some detail and the elements of the asymptotic theory were worked out. 
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A subsequent work of Lees ( 1947 ) carried out calculations of neutral- 
stability curves up to = 1.3, and also included the famous prediction 
that cooling the wall would' act to stabilize the boundary layer* The 
next study was by Dunn ( 1953 ) and Dunn and Lin (1955 )> where the asymp- 
totic theory was worked out in a less approximate form than by Lees and 
Lin, and three-dimensional disturbances were also considered. Neutral- 
stability curves were computed at = 1 . 3 , 1*6 and 2,2 by Mack (i960) 
on the basis of the Dunn-Lin theory, with the inviscid solutions being 
obtained with improved accuracy, The next step was by Resliotko (i960) 
and Lees and Reshotko (1962). The asymptotic theory was retained, but 
both the inviscid and viscous solutions were obtained by numerical inte- 
gration. The viscous equations used contained higher-order terms than 
in the usual approximation. 

Direct numerical integration of the compressible stability equations 
was first carried out by Brown (1961, 1962), first for the simplified 
equations of Dunn and Lin and then for the complete linearized equations. 
More recently, Brown (1967) has integrated a more complicated system of 
equations which include the mean normal boundary -layer velocity and are 
valid for three-dimensional disturbances, I started my own work on the 
numerical integration of the stability equations in 1961. The numerical 
method is described in Mack (1965a ), and results for two-dimensional dis- 
turbances are presented in Mack (1965b). 

13*2 Linearized., parallel flow equations 

A comprehensive account of the theory must start, as it did in Part 
A, with the derivation of the governing equations. The dimensionless 
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Navier- Stokes equations for a viscous, heat-conducting perfept gas 'are 
as follows . 



A 


■3kj 





{13'. i) 



(13.2) 

(13.3) 




(13A) 


where 



( 13 - 5 ) 


The summation notation has been adopted for economy, • and the asterisks 
refer to dimensional quantities- The equations are, respectively, 
momentum, continuity, energy and state. The quantities which did not 

w. n 1 * 

appear in the incompressible equations are T , the temperature; ,k , 
the coefficient of thermal conductivity; R , the gas constant; ct , 
the specific heat at constant volume, which will be assumed constant; 
and X , the coefficient of second viscosity {= 1.5* bulk viscosity 
coefficient). 

The stability equations are obtained from the Navier-Stokes equations 
by the same procedure that was used for incompressible flow. All quan- 
tities are divided into mean flow and disturbance terms, the equations 


- 98 - 


Ut.GEMAL PAGE JS 
'OF 



are linearized, the mean-flow terms are subtracted out, and, finally, 
the parallel -flow assumption is made* The resulting equations are 

i 

14. 

then made dimensionless with respect to , the velocity sca3e, S , 
the length scale, and free -stream quantities serving as reference 
quantities for all state variables and transport coefficients. ! '- The 

* t H 

exception is £ , which is referred to Cp^u, . The x-momentum equation 
is ^ 



^ t AT* Aj, J (13. 6) 


The y-monientum equation is 
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( 13 . 7 ) 



The a -momentum equation is 
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The continuity equation is 
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The energy equation is 
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The equation of state is . f 

f'r ^T' ■+ o*T • (I3.ll) 

Additional quantities which appear in these equations are Mi , the 
free-stream Mach number; V , the ratio of specific heats; and v ** c ? ^ 
is the Prandtl number and is assumed constant- • 

The boundary conditions at aj. = 0 are 

= c y v't*) a 0 f w'Cal «. 0 , T'c*) - 0 ' (13.12) 

The condition on the temperature fluctuation is suitable for a gas 
flowing over ‘a solid wall. For almost any frequency of the gas, it 
is 't possible for the wall to do other than remain at its mean 
temperature. The conditions at infinity are that as -e» 00 , the 
disturbances must at least be bounded.- This condition is less restrictive 
than requiring that the disturbances go to zero. In supersonic flow, 
waves may propagate to infinity and we wish to include those that do 
so with constant amplitude. 

13 . 3 Reduction of equations to near two-dimensional form 
We now specialize the disturbance to a wave of the same- form as in 
Part A. 

e*p[i * + gi 4 - tfs C £)] (13.13) 

With the disturbance in this form, we inspect (l3.6)-(l3.1l) to see to 

1 

what extent they simplify in the tilde coordinate system, i.e*, the 
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coordinate system where the z axis is parallel to the wave front. It 

■w *-r * ^ 

is evident that with 3/33 -o , all w and W terms drop out of the x 
and y momentum equations and the continuity equation. When, further, 
the reference velocity is changed to U ,* cos \y , where V" the angle 

*i ' * 

between the wave normal and li t 3 these three equations are identical 
to the equations for a two-dimensional disturbance in the mean .flow 

with Mach number Mj - Mi co$ and Reynolds number c°s . 

If the boundary layer in the original x, 2 coordinates is two-dimensional, 
then However, there is an important difference here from 

the situation in incompressible flow. In the latter case, the transformed 
flow corresponded to a real flow. Here, is the boundary-layer 

velocity profile for the free-stream Mach number not 1*L . Consequently, 
even if the energy equation can be reduced to two-dimensional form, the 
tilde equations do not represent a real flow. The tilde coordinate system 
thus serves only to simplify the equations, not, as in the incompressible 
case, to relate all three-dimensional disturbances to the solution of the 
same flow for two-dimensional disturbances. 

It is evident from an inspection of the energy equation, ( 13.10), 
that even with 9/3?= 0 all terms which include v’and W do not drop out, 

*V“ 

Two dissipation terms with w and W remain. This result, that the three- 
dimensional stability equation cannot be transformed into two-dimensional 
equations, has been known to most workers in stability theory for a long 
time although it is rarely mentioned in the literature. I first learned 
of it from Prof. M. Finston at M.I.T. in about 19 ^ 8 . It is mentioned by 
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Dunn and Lin (1955) in the introduction of their paper., but the first 

explicit demonstration seems to be by Reshotko ( 1962 ), This trans- 

* 

forma' ion is usually discussed by working with the equations for the 
amplitude functions and the physical significance is obscured by the 
manner in which the derivation is carried out* The type of derivation 
used here was first given by Dunn and Lin (l955)j although they too 
reverted to the older derivation to obtain their final equations. The 
complete derivation along these lines was first carried through by Shen 
(1964 )j but only for the simplified equations of Dunn and Lin where the 
dissipation terms of the energy equation are neglected. 

The dissipation terms of the energy equation are 




Alt. f 9ti # 


+ tt 


Xu rain 1 . a* 


<tr 


T , + 1« d* •£' 


The second group of terms has a term with the mean velocity "W , and this 
term introduces no essential change from the two-dimensional equations. 
For a two-dimensional boundary layer, where - UO) , the second 

group of terms can be written 


in r (*if „ A*) 1 ] T < = di (&X II 

tlT \Xifl l"* 1 }'. <<T\*t^] c»s l y (13.15) 

However, nothing can be done' with the final term, which involves w' , 
the fluctuation parallel to the wave front. That such a fluctuation will 
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exist is evident from the a -momentum equation. 

t * * „ (13.16) 

. + A^ rr AV + 1$L + ?!' 4. tL^T ; ) 

at 1 ^ ^ *V * 

We cannot have $?'= a as a solution, partly because of the source term 
v f A*W / Xt^ . There is always a source term of this type whenever a 
normal-velocity fluctuation, v' , is combined with a normal gradient 
of a mean-flow quantity. An example in the Orr-Soramerfeld equation 
is v' A U./ ({j* , which is a vorticity source term ( Ah( is the mean 


vorticity). 

With the i?' term included in the energy equation, it is necessary 
to use the 2 -momentum equation, and, as we shall see, the order of the 
system, which is six fo* a two-dimensional disturbance, becomes eight, 
Dunn and Lin (1955) used a simplified system of equations in which all 
of the dissipation terms are neglected. When this is done, the eigen- 
values of three-dimensional disturbances can be found from the two- 
dimensional equations. However, it is still necessary to use the z- 
momentum equation if the eigenfunctions of the velocity fluctuation 
components in the x-z plane are to he calculated. An important example 
of equations with no dissipation terms are the inviscid equations, 

13.4 Introduction of Fourier components 

With the fluctuations in the form (l3»13)> the complex amplitude 
functions for u r , v' , ^ and T' are defined to be ? , * ^ , 

X 3 T , and 8 . (Since v' * v'/Wf and v' a & cj> ; ^ cos ^ ■ , 

we have &' = 3 <3 )* The equations for the amplitude functions are 


* 
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pluK-cU * -Sr'j] = 
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where 


s s 

JiT 


(13,23) 


The subscript S has been dropped from and f? . If, whatever length . 
scale, C , is used, ( )' is interpreted to mean differentiation with 
respect to ^ / L* , then the corresponding wave number and Reynolds 
number become 05 } R u * * In particular, if 



X 4 


where R* is the x-Reynolds number 




j then 



R = ik _*L 

M* (fit ) 1 ' 2 



1 



(13. 2^) 


(13.25) 


( 13 - 26 ) 


Consequently, with the y variable replaced by the Blasius similarity 
variable vj , then k is given by (13*25) aud R is simply . 

This system permits boundary-layer profiles calculated in terms of the 
similarity variable to be used directly in the stability equations. 

However, the boundary-layer thickness, displacement thickness, or momentum 
thickness can be used equally well as the reference length, with the proper 
interpretation of , R , and the independent variable. All of these 
<&$ and Rs differ only by numerical factors which depend on the Mach 
number (e.g., ^ ~ 7* od , = J h + cl , R# = >?<> R )• Since the numerical 
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results are to be presented in terms of <* and & , Table 13 »1 is included 
to give these numerical factors for the instated flat -plate boundary layer* 


% 

Vs 

V... . 

V* 

0 

6.0 

1.72 

0.664 

0.7 

6.2 

1*92 

0.660 

1*0 

6.4 

2.13 

0.656 

1. 6 

7.0 

2*77 

0.648 

2*0 

7-6 

3*37 

0.644 

2.2 

8.0 

3.72 

0.643 

, 3.0 

9*8 

5*48 

0.642 

3.8 

12*1 

7.83 

0 . 644 

4.2 

13*5 

9.22 

0.645 

4*5 

14.6 

10.34 

0. 646 

4.8 

’ 15.8 

11.55 

0.646 

5.8 

20.0 

15.73 

0.636 

6*2 

21.7 

17.49 

0,629 

7.0 

25.4 

21.19 

0, 616 ' 

7-5 

27*8 

23.62 

0.607 

8.0 

30.3 

26.13 

0.598 

8.5 

32*9 

28.72 

0.590 

9.0 

, 35*5 

31.38 

0.581 

9-5 

38.2 

34.10 

0*573 

10.0 

4i,o 

36.88 

0.565 

13.1 

Boundary-layer thickness (u x / li? = 0.999),? 

displacement 


thickness and momentum thickness of insulated-wall boundary 
layer. (Wind-tunnel temperature conditions) 
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13.5 System of first; -order equations 

In order to carry out .the numerical integration of the compressible 
stability equations , it is convenient to reduce them to a system of first- 
order equations. Following Lees and Lin (1946)., we define the dependent 
variables for this purpose as 



(13.27) 


The fact that this choice of dependent variable leads to eight first-order 
equations demonstrates that we have an eighth -order system. The £ 3 

equation is obtained from the continuity eouation; the equation 

***■•)■ , 

comes from the x -momentum equation; the £ 4 . equation, after a lengthy 
manipulation from the y -momentum equation; the £5 equation from the 
energy equation; and x^ equation from the z -momentum equation. The 
resulting eight equations can be written as 'the product, of a matrix 
( d;j ) arid a column vector £ . 





(13.28) 


t - 1 


The 30 non-zero elements of the 8 x 8 matrix ( ) are listed below, 

n 1 

The equation has only one non-zero coefficient. 



(13-29) 
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l 

In the z l equation, six coefficients are non-zero. 




~ i d. R It - c - z 
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(13.30T) 

The 2 3 equation has four non-zero coefficients 

' ' 
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The 2 a equation has six non-zero coefficients. The coefficients 
written in terms of the factor L , defined by 
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The Zj- equation has only a single non-zero coefficient, 
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(13.33a) 

(13.33b) 

(13.33c) 

(13.33d) 
(I3.33e) 
(13-33 t) 

(13.3^) 


-no- 


.1 


is/ i 


The equation has six non-zero coefficients. 
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The 2, equation has only a single non-zero coefficient* 
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The 2 y equation has five non-zero coefficients. 
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In these equations, 

4=A 

/“ 

and is a constant* 

The hound ary conditions are 


(13.38) 


2 , (»' = 0 t = 0 , Sj-f 0 ) » 0 f 2 , Co > = 


(13-39) 


and at ^ , the amplitude functions are hounded. Just as for incom- 
pressible flow, we have to deal with an eigenvalue problem and can only 
satisfy the boundary conditions for certain combinations of «t , | , c 
and R . 


13*6 First-order equations in free stream 

In the fre^ -stream, yj > , the equations reduce to the following 

equations with constant coefficients. 


t 

2, = 2; 


( 13 . 40 ) 
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(13. Va) 
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The general solution of these equations is of the form 
3 


*1 - E c i A'i'e^-y" , 


I a I, ... , 8 


i = « 


(13.45) 


(13-46) 


(13.47) 


The Xj are the eight characteristic values of the coefficient matrix, 

C I the A** are the 64 components -of the eight characteristic vectors 
A ^ 1 ; ana the eight C ^ 1 are arbitrary constants to be determined by 
the eight boundary conditions* Equations (l3*40)-(l3.46) are simple 
enough so that the Xj and A ^ 1 can be obtained analytically. 

A derivation of the Aj and A; is given by feck ( 1965 ) for the 
sixth-order system. The extension to the eighth-order system is obvious. 
Of the four characteristic values which satisfy the boundary conditions as 
17 -r to t only two have simple forms. They are 

.Vi 


A - -*[d + I cl R O - cA| 


(13.48) 


- 113 - 



and are identical to A 3 of (4.17), the characteristic value .of the 
viscous solution for the Orr-Soramerfeld equation in -the free stream. • f 
In the limit of large &R the four characteristic values are, vith the 
numbering chosen to correspond to the incompressible cas$, 

X, S - ol [l - H - cVj * (13.49) 
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To the same approximation, the characteristic vector A has the com- 
ponents (with Ai 1 , A ( f and A^ 1 omitted in the listings) 
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The characteristic vector has the components 
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The characteristic vector A has the components 
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(13-554) 




Finally, the components of the characteristic vector A are 


A?' = A?' - fu 


IM 

A* 


(13.56a) 


A, = 1 


(13.56b) 


The first of these vectors is the compressible wavy -wall inviscid 
solution; the second is the viscous velocity solution; the third is the 


- 115 - 


viscous temperature solution; the fourth is a second viscoup velocity 
solution for the non-periodic (in 5 ) velocity component parallel to 
the wave front. The exact counterparts of these solutions provide the 
initial values for a numerical integration which proceeds from the free 
stream to the wall. 

A discussion of the inviscid equations in the free stream will he 
given at the end of the next section, 

13.7 Inviscid equations 

When all of the terms with viscosity and heat conduction are dropped 
from ( 13 . I?) “(13 -22), the equations reduce to 


f [icff-o ? + £}] ^ — 

ST M, 2 - 

(13.57) 
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p [i ( U- <■)(> * f'j] = t fj 

(13.61) 

ir = a t i- 
? T 

(13.62) 


In these equations,, as already indicated in Section 13*3, and 
W appear only in the z-momentum equation, and the other equations are 


identical to two-dimensional equations* Consequently , we cam obtain 
the eigenvalues of all inviscid three-dimensional waves frpm the two- 
dimensional equations, although in contrast to incompressible flow, 
separate eigenvalues must be obtained for each wave angle , The 

amplitude function can be computed from (13*59) after the eigenvalues 
are known without an additional integration. 

13 . 7*1 Various forms of inviscid equations 

We can easily derive several forms of the inviscid equations which 
will be useful for various purposes. Substitution of (13.62) into (13.60) 
yields, with the use of (13.61), 

l V ( i i + <P ' ) 

IT = ±~-t 

li- c 

When this equation is differentiated and the result substituted into 

* 

(13*58), and *f eliminated from (13-57) and (13*63), the result is 
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( 13 . 64 ) 


This is the form of the invscia stability equation derived by Lees and 
Lin ( 1946 ). Since it is a second -order equation, just as for incom- 
pressible flow, only two boundary conditions may be satisfied. 
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Equation (13.64) is useful for analytic developments, but it (is 
not suitable for numerical integration because of the denominator 

(w l w 1 

T- , which, for C; so , is zero at the point in the boundary 

* 

layer where U- is equal to the local speed of sound. From (13-57) • 


and (13-63), 

5 . T - S' (it - c) 

* " ‘ t- (13.66) 
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Substitution of this equation into ( 13 . 63 ) gives 
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( 13 - 67 ) 


When ( 13 . 67 ) is solved for flj>’ , we find, with the 2 notation. 
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(13*68) 


and the y-momentum equation (I3.5S) is 



(13.69) 


These equations are the counterparts of (5.14) and (5*15) in Part A, 

When M, = 0 and T= 1 f they reduce to the incompressible equation. 

For some purposes it is useful to think of an observer moving with 
the wave. In order to write the equations in this coordinate system, we 
introduce 
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For Cj ■= a j M r is the local Mach number of the relative flow. Then 
(13.67) becomes , 


Tr 
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(13.71) 


This equation can be recognised as the linearized pres sure -area relation 
of one-dimensional flow. The quantity cp / (1C- c) is the amplitude function 
of the streamline slope. 

The other disturbance amplitude functions can also be written in terms 
of Q l . However, because of the mean normal gradients, the complete 

expression for the amplitude functions will also involve source terms. 

We find 


r~- 

f 



li-c 


1 , 

CTi, { J/C — C- *1 

1“ Mr 



(13.72) 
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(13.73) 
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1 - m‘ ^ 

k- J. 


(13-7^) 


When the second terms of these equations are written in terms of Tr from 
(13,71), they can be readily recognized as the linearized momentum equation, 
isentropic temperature-pressure solution, and isentropic density -pi cssure 
solution, respectively. 
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We are able to derive a single second-order equation for ^/(Zt-6) . 

It is ’ ' 

(13.75) 

A similar equation can be found for TT . 



Tr" - (l n Mr Tr* S l (i- M^TT <y 

u 


(13.76) 


This equation has been studied by Lighthill (1953) in another connection. 
13.7*2 Inviscid equations in free stream 
In the free stream., (lj. 76 ) reduces to 

1 r' - **0" “ 0 (13.77) 

•v ' 

where Mr, is a constant and is given by 

M Ri = Ci-cl M, (13.78) 

The solution of .(13*77) which is bounded at infinity is 

Tr = - * M r , cvp [- 3 (l- iy - ( 13 - 79 ) 

vl !/l 

The characteristic value, - ( 1 - f is equal to the X, of Section 

13.6. In order to simplify the notation in the remainder of this section, 
Mr ( will be replaced by M a and ct by ct . 


- 120 - 



When C; s o and M a. A I j *ir o as y , We pall this ease 

a subsonic neutral disturbance. When M ^ > l } we have a supersonic 
neutral disturbance. The coefficient of in (13*79) is pure 

imaginary and the wave amplitude is independent of y * The whole 
solution is 


lmI- ir<7'7im (13.80) 


It is evident that these waves are Mach waves. There are two cases 

l /r * 

depending upon the branch of the square root we choose. If (-0 - t I ; 

we have the minus sign in (13.80) and the case shown in (a) of Fig. 13.I. 
If (-0 1 =■ - i j we have the plus sign and the case shown in (b). 



(a) (b) 

Fig, 13.1 Supersonic neutral disturbances 


Case (a) is an outgoing wave; case (b) is an incoming wave. Case (a) 
appears more appropriate to the bound ary- layer stability problem. However, 
a combination of the two types of waves permits the boundary conditions to 
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be satisfied for any * and c and ve no longer have an eigenvalue 

* i 

problem. But if we restrict ourselves to only one kind oj* wave, then 
we still have an eigenvalue problem. Type (a) waves havQ actually been 
found, as will be seen in Section 15-6.3* 

For amplified and damped waves, c; £ o and Mr will be complex. 

M r r O-tlMi » (i* O Mj 

( 13 . 81 ) 


The res! part of Mr is always positive, but the imaginary part will 
change sign depending upon whether the wave is amplified ( c ; > 0 ) or 
damped (c; < 0 ) . We can write 


z \ L r ■/ 



1- (MrC. 


(13*62) 


under the assumption that (Mr); *c < (Mr)* . When (M^ h < 1 , the 

boundary condition at infinity is always satisfied. The wave fronts of 

the amplified disturbances are tilted in the upstream direction; the 

wave fronts of the damped disturbances are tilted in the downstream. 

direction. (The neutral wave fronts are normal to the free-stream 

direction.) When (M r )|. > 0 , there are again two cases to consider. 

*( L 

With f'-i) - + 1 , we have the outgoing case (a) waves: 


exp [** c< (1 - Mr) \7~7sl] 


oi (Mr]* (MrI I 
I — — r~y-77r- li 


• [(Mr)* ~ l] 


'fl 


, f X 'll . 

[(M N )r j] (13.63) 


-122- 



v* . . 

With (~i\ = -i , we have the incoming case (b ) waves. 


exp [-*(*- Mr] (7-7^] 


l 


c< fMnV (Mr); 


=• e *P S- — 1 {'h - 7 ^ 


t(M„C - 0 


<2X 




(13-84) 


We can classify these waves by whether they increase or decrease with 
increasing 7 . 

Amplified Damped 

Case (a), outgoing ’ J decrease Increase 

Case ‘(b), incoming increase decrease j 

For amplified waves, only the outgoing waves satisfy the boundary 
condition at infinity. For damped waves ( C\ ^ o ) } only the incoming 

waves satisfy the boundary conditions. Eence only two of these solutions 


can be accepted. 



l4. Numerical Integration of Stability Equations 


14.1 Inviscid equations 

Two methods have been devised for the numerical integration of the 
inviscid stability equations. The first of these is by- Kqshotko (i960). 

The second-order linear equation is transformed into a first -ordex' non- 
linear equation of the Riccati type. This equation is solvec by numerica. 
integration except for the region around the critical point, where gener- 
alisations of Tollmien's series solutions are used. Although the numerical 
integration is performed on a computer, the procedure is not completely 
automatic because of the matching required between the numerical inte- 
gration and the series solutions. 

The second method is by Mack (1965a) and is a. generalisation to com- 
pressible flow of Zaat’s (195S) method. The numerical integration is 
carried out along an indented contour in the complex plane. The compu- 
tation of the eigenvalues is completely automatic provided the initial 
guess is adequate, a condition that can easily be met. The equations 
that one uses are the two first-order equations (13.68) and (13. 69). The 
integration starts in the free stream and continues to the wall. In the 
free stream, the equations can be written as a second -order equation for 
2* as has already been done in Section 13*7* In the fixed coordinate 
system this equation is 





[\ - KO - cf] a 0 


(l4.l) 
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The solution which satisfies the boundary condition at infinity is 


= - Cl- cl exf 3 [l - m‘ C.I- O*] 


(14.2) 


The integration constant has been chosen to make t ^ 51 C *, o ^ . 

The initial conditions at '‘j - for the numerical integration are 


ity = i;h~ 


>V'« 


(14.3) 


2 4 = - n- a 


(14.4) 


We must decide which of the indented contours shown in Fig. 5*1 
to use. For a neutral disturbance., both contours are satisfactory. 

For a damped or amplified disturbance, the choice of the correct contour 
is made on the same basis a$ in incompressible flow for a velocity profile 
with an inflection point (p. 3*0* Table 1^.1 gives some eigenvalues for 
a boundary layer on an insulated flat plate at = 3*8. Contours (a) 
and (b) are those of Fig. 5*1* 


Contour 

0 ! 

Cv 

c i 

a 

O.I787 

0.8585 

-0.005 

b 

0.1053 

0.8317 

- 0.005 

a 

0.1058 

0.8317 

+0.005 

b 

O.I787 

O.8685 

+0.005 


Table lt.l Inviscid eigenvalues at M_ = 3.8 for different 
indented contours ■ 
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The eigenvalues of the neutral disturbance are oi 5 = 0 . 1421 , - O.8523. 

At finite Reynolds numbers, the unstable region lies below the neutral- 
stability curve; the stable region above it. Consequently, the criterion 
that the inviscid solutions must be the R (=0 limit of corresponding 
viscous solutions selects contour (a) just as it did foy incompressible 
flow. 


The numerical integration yields for specified values of o[ , 

Cy. and c; . To find the eigenvalues, one of the three parameters is held 
constant and the other two perturbed until the boundary condition M - 0 
is satisfied. With a reasonable initial guess, a linear perturbation is 
adequate. If a neutral solution is sought, then c; r 0 , and ot ■ and c K 
are perturbed. An initial integration, followed by two more integrations 
with ot then cv perturbed, yields the following four quantities;. 


A 2^ M 


A f 3; M 


Act 


A ct 


A £3^ Co\ 


AC, 


A ( 0 ) 
ACv 


The new t X and Cy. are found from 


— — A ct -t- A Cy 




A* 


AC, 


^ ^3,-^ ^ c< + & to} A c r „ 2^ ( 0 } 

At* A Cy 


( 14 . 5 ) 


( 14 . 6 ) 


( 14 . 7 ) 
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This process is repeated until and c y converge to the specified 
accuracy. If oi is held constant, then one integration cab be saved 
by taking advantage of the fact that Zj and Z A are analytic functions 
of c , The two derivatives similar to those of (l4.5) tfhich are not 
actually computed are obtained from the Cauchy-Riemann equations. Since 
there is no serious problem of error buildup with the inviscid equations, 
single-precision arithmetic and a standard forward-integration method 
are adequate. 

Ik . 2 Complete equations 

The numerical integration of the complete compressible stability 
equations presents the same central difficulty as with the Orr- Sommerf eld 
equation. Two methods have been successfully applied to this problem. 

Both methods were first restricted to two-dimensional disturbances, but 
have since been extended to also handle three-dimensional disturbances. 
Double-precision arithmetic is used in both methods to control the error 
build-up. The first method to be developed, just as was the case for the 
incompressible boundary layer, is by Brown {1961, 1962). In Brown's 
method, three arbitrary independent solutions (for the sixth-order 
system) are produced by numerical integration -from vj = 0 to 
At ^ , these solutions are matched to the analytic free-stream- solutions. 
Two of the matching conditions can be immediately satisfied, but the third 
condition can only be satisfied by searching for a suitable combination 
of cl , c and R . Brown successfully applied his method to several 
problems. 

The second method, which was developed independently shortly after 
Brown's method is by Mack (1965a). This method is similar in concept to 
the method described in Section 14.1 for the inviscid equations. The 
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analytic free -stream solutions provide the initial conditions for four 
linearly independent solutions (eighth-order system). These solutions 
are produced by performing four numerical integrations from ^ to 
•vj = o . Three of the four boundary conditions at ^ - o can be im- 
mediately satisfied. The remaining condition, (o\ .o , is satisfied 

by the same sort of linear perturbation procedure described in the pre- 
ceding section. With ci and R held constant, only a single pertur- 
bation integration (of each of the four independent solutions) is required 
because the it are analytic functions of c .To calculate a neutral 
solution Ccjso'i jC*. is held constant and and R are perturbed. 

As is often the case in problems of this kind, . the fact that H 5 - C o \ is 

a physical quantity leads to a rapid convergence. Often only a single 
iteration is needed to give the eigenvalues to three significant figures, 
and it is rare that more than two iterations are required. Once, the eigen- 
values are known, then a subsequent integration of a single solution which 
is formed from the four independent solutions : in the proportion 'determined 
during the eigenvalue computation will produce the eigenfunctions and 
whatever related quantities are desired.. 

On an IBM 709 ^ computer, this method requires 0.25 sec to integrate, 
the four solutions of the eighth-order system across one integration step; 
O.lk sec are required for the three solutions of the sixth-order' system. 

With 160 steps, a common number, and provided only one iteration is needed, 
it takes tx^o minutes to calculate the eigenvalues of the eighth-order system, 
and 67 sec. for the sixth-order system. At a computer charge of $250/hr., 
the costs are, respectively, $8.00 and $UA5 for each set of eigenvalues. 
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These numbers are included here to emphasize that the economics of doing 
research on a computer is a major consideration. 

The method just described has one important drawback* It is limited 

t 

to maximum values of otf\ which are not always adequate. Fortunately; as 
increases, so does (tf. R , but a value of 300 is about the highest one 
can achieve. The adoption of some version of Kaplan's method; such as 
the method of Wazzan, Okamura and Smith (i960), would remove this limitation. 
Also since single -precis ion arithmetic could then be used, the speed of the 
program would also increase. 



15* Inviscid Theory 


15*1 General results 

Lees and Lin (1946) have carried out the most complete study to 
date of the inviscid theory, Their results which are useful for our 
purposes will he summarized here without proof. Since the equations 
for a three-dimensional disturbance can be reduced to the two-dimensional 
equations, only the latter equations will be considered in this section. 
Lees and Lin classified the disturbances into three groups, depending 
on whether the phase velocity c£ is subsonic, sonic or supersonic 
with respect to the free-stream velocity U? . In dimensionless terms, 
this classification is 


C*. > i - 


Mi 


subsonic 


Ci. = l- 


Mi 


sonic 


(15.1) 


c, < 1 - ir 

Mi 


supersonic 


Lees and Lin's chief results are;' ' 

(i) The necessary and sufficient condition for the existence of a 
neutral subsonic disturbance is that there is some point in the 

boundary layer where 

A f* 


A U-\ 

Ay It i 0 


(15.2) 


The phase velocity of this disturbance is c s , the mean velocity at 

/ t a. 

. (The proof of sufficiency depends upon T- Mj tk-c) ><?.) This 
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condition is the generalization of the condition found in Section $.1 
for incompressible flow that there must he an inflection point in the 
velocity profile for a neutral disturbance to exist* The ppint y s , 
which plays the same role here as the inflection point, will be called 
the generalized inflection point. 

(ii) A sufficient condition for the existence of an amplified 
disturbance is that , 



at some ^ greater than , where is the point at which C- i - */m,. 
The proof of this condition also depends upon T~ M? { U~ <c*l > 0 , 

(iii) A sonic neutral disturbance is t 

ot te o , Co - 1 - — - • (1J5.3) 

M, 

(iv) If T - M, 1 (H- c'i > 0 t there is a unique wave number e{ s 

corresponding to c s for the neutral subsonic disturbance. In other 
words, as long as there is no region in the boundary layer where the 
local Mach number of the mean flow relative to the wave velocity is 
supersonic < 1 ) . , a unique exists. 

These results are obtained by a direct extension of the methods of. 
proof used for incompressible flow. The necessary condition for a 
neutral subsonic disturbance is derived from the equation for the 
discontinuity of the Reynolds stress - ^ uV' at the critical point. 
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This equation is, for Cj = o , 

|e? f 

- tCy s -o) = 4 «T ilii 

1 1 k: 

Since % is zero at the vail and in the free stream by the boundary 
conditions for a subsonic disturbance, it follows that Cu'/r) must 
be zero at . A quantity that appears in the asymptotic viscous theory 
is 


14 /V 

,V T 


J 7‘ 


(15.4) 


= — 7r 1 (hj 

Tv, (u:f 3 [V T 'J, 


(15 -. 5 ) 


In terms of v<.CcY , 


„ i Tw [lQc lli'f 


(15.6} 


The discontinuity in 'f is derived from the generalization to 
compressible flow of the Tollmien series solutions for (J>, and 
These solutions are 




(15.7) 


^£•7! = jity >' 


T, 


s r 


Od\\*j\rn 


e! (U‘ 


1^(7-^) , (15.8) 
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For iy < s - i tt as for incompressible 

flow. The leading terms of and are U't. and X / lit , re- 
spectively, so that and Gj t are normalized here in a different manner 
than in Section 5*2.1. These solutions have been worked out in more 
detail by Reehotko (i960). Both tp and $ have the same analytic behavior 
as for incompressible flow. What is new here is the existence of a temp- 
erature fluctuation. According to Reshotko, it has the behavior 


X 


4 (u‘ 


-_jj Uj.( 7 -v ,«1 


•I- * * * 


nil ,M 

* 7 C 

Hence, even for a neutral disturbance, 'here fy! / T) ~ 0 and ^ and 
■f are both regular, 6 has a singularity at *^ c . 

The inviscid solutions in powers of ot 2, can also be extended to 
compressible flow. They are, according to Lees and Lin (19^6)^ 


(15*9) 


03 

g, . W til- c\ J2 


T»a 0 


( 15 * 10 ) 


O0 

= Ck-c)^ k^c-/) 


o 


where 


M 0 {y\ & 1 

-5 

« [ 


T- M Uu~c) 1 f f ClL-c) 


(i9.ll) 


Oi-d 1 
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and 



. / "T*~ Mi f O 
«4 Ot-c f ' 



T- MhCZt-cf 
_ Ot-c? 




( 15 - 12 ) 


In the course of my numerical studies of the eigenvalues of the 
complete stability equations,, it became evident that (jVl could not 
be true for T- Mi Cli-c) ^ o ( ^ i) , The consequence of this loss 

of uniqueness makes the inviscid theory for supersonic free-stream Mach 
numbers differ greatly from the incompressible theory® 

Before going further, we must examine more closely the consequences 
of the finding that neutral and amplified inviscid disturbances can exist 
whenever (jJ- / T ) ~ 0 „ For the flat plate incompressible boundary 

. |t 

layer. It is negative everywhere except at y ~ o „ However, for a 
compressible boundary layer on an insulated flat plate (it/T)* is 
always zero somewhere in the boundary layer. Consequently, all such 
boundary layers are unstable to inviscid disturbances. We have seen in 
Section 9-2 for the Falkner-Skan profiles with ^ < o f that inviscid 
instability became more important as the inflection point moves away from 
the wall. Figure 15.1 shows <= 3 f the mean velocity at the generalized 
’ inflection point and the phase velocity of the neutral subsonic disturbance, * 
as a function of M-^, The wall is insulated, and the free-stream temperature, 
which is a parameter for the exact numerical solutions of the boundary-layer 
equations used here, is characteristic of wind tunnel conditions* The 
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1.0 



Fig» 15.1 Phase velocity of neutral subsonic disturbance 
(upper curve) and of neutral sonic disturbance 
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* o 

stagnation temperature is held constant at 311 K until, with increasing 

Mp the free-stream temperature Tj* drops to 50°K» For higher Mach 
numbers, T k a is held constant at 50°K- All results to be presented for 
insulated -wall boundary layers are for this family of boundary-layer 
profiles unless otherwise noted. 

Figure 15*1 also includes a plot of c 0 s= 1 - i/m, « For a 
disturbance to be subsonic relative to the free stream, and hence have 
vanishing amplitude at y even for a neutral disturbance, must 

be greater than c„ , It is often said that only subsonic disturbances 
are consideied in stability theory, a statement that is not entirely correct 
It is true that the neutral subsonic disturbance ot 5 only exists when 
C £ > \ - ‘/M.* . However, this doesn't rule out amplified or 

damped disturbances with < I " '/M, , or even neutral supersonic 

disturbances with a c h different from. c 5 , We shall find examples of all 
of these disturbances, all of which satisfy the boundary conditions at 
infinity and lead to eigenvalue problems. For c* y- 0 } the amplitudes 
of outgoing amplified and incoming damped waves vanish at infinity re- 
gardless of the value of • for cj » 0 , the amplitude wi31 only be 

bounded at infinity for < c„ . What does turn out to be true is t’ at 

the most unstable disturbances are always subsonic. Further, we shall find 
that for one class of disturbances, the amplified first-mode waves, 

C 0 ^ <■ Cj , This result has important consequences, 

15*2 Multiple neutral solutions, two-dimensional disturbances 
15,2.1 Neutral subsonic solutions 

An important omission in the Lees -Lin inviscid theory was the failure 
to recognise that for a boundary layer at supersonic free -stream Mach 
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numbers , there are multiple neutral solutions. It is rather easy to 

, '* 1 .- 

aemonstrate the existence of these solutions. The inviscid equation 

■written in terms of tt , (13*76)* quite evidently has a different 

character depending upon whether is less than or greater than unity. 

It is instructive to consider large enough so that the TT* tern can 

be i glected. Then ( 13 . 76 ) reduces to 

rr“ - M^)7T - o (15.13) 

When <. [ , the solutions of (15.I3) are elliptic in nature, and 

it ±f under this circumstance that Lees and Lin proved the uniqueness 
of . However, when > t , (15.13) becomes a wave equation, and 
as in all problems governed by the wave equation we can expect there to 
be an infinite sequence of wave lengths which will satisfy the boundary 
conditions. 

If we introduce 


v = ( V\ 


1 


then , from ( 15 . 13 ) * 


T = ! cos 


A X Ay 


( 15 * 14 ) 


7 <■ -j, (15.15) 


TT 




7 >■ (15.16) 


where (15*15) follows from (13.69) and the boundary condition on at 
1 

Wot too long after I found these solutions for the boundary layer, 
Gill (1965 ) independently found them during Ms study cf the "top-hat" 
wake. 
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^ = O ' . We have written oc as . The subscript 5 designates the 
neutral subsonic solution; the subscript h refers to the multiplicity 
of solutions. Since we have arbitrarily made T positive for y ■> , 

it can have either sign at j a a # At } * J, ; 


COS 


r * 

^5-n / ( - i) Ay 


=i t 1 


(15 „ 17 ) 


and 



h s l, ij . . ■ (15.18) 


Because of the approximate nature of (l5-15)j the magnitude of tf 5Vl as 
given by (I5.l8) is of little importance. However, the difference between 
adjacent values of d ST / , given "by 


. (15.19) 

tMft- t\' 1 »t 7 

turns out to be accurate under many circumstances. 

When the numerical integration of (13*68) and (13.69) is candied out 
with ss c s , C; * (f for the insulated -wall boundary layer, the 
which are found by the eigenvalue search procedure described in 
Section li*-.l are shown in Fig. 15*2. The solution for each value of n 
will be referred to as a mode, with vi — 1 the first mode, h -1 . the 
second mode, etc. The wave numbers of the first mode were first computed 
by Reshotko (i960). The second and higher modes appear only for > 2. 

t 

With £V = c a and the family of insulated -wall boundary-layer profiles, a 
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point *j ( at which = 1 first occurs at - 2,2 where y, f 0. With 
increasing the point moves out into the boundary layer- No higher 
modes could he found numerically for 2.2, in accord with the theory 
given above. Also, from (15.18) c< 4h must be inversely prpportional to 
ij, , the thickness of the relative supersonic layer. As M^-^2,2 from 
above, 0 and must become infinite, again in agreement with the 
numerical results- ' 

A curious feature of Fig. 15.2 is that the upward-sloping portion 
of the first -mode curve between « 2 and 4,5 is in a sense continuous 
through the other modes; i.e., there is a Mach number range for each mode 
where the vs M-j_ curve has a positive slope. The end point of this 
region for one mode is close to the starting point of a similar region 
for the next higher mode. The approach becomes closer as increases. 

Thio upward-sloping portion of each curve has a special significance which 
will be pointed out later. The same phenomenon will appear in even more 
emphatic fashion when we take up three-dimensional disturbances in Section 

15. 5. 

The accuracy of the wave number spacing given by (15.19} can be checked 
from the numerical results, which are given in Table 15.1 for =5*8 and* 
10 . 0 , 

« 5.8 U x * 10.0 


Modes 

a oc 


2-1 

o.i4 

0.09 

3-2 

0.13 

0.10 

4-3 

0.27 

0.08 

5-4 

0.28 

0.02 

6-5 


0.10 

7-6 


L 10 

8-7 


0.10 

(15.19) 

0.28 

0.10 


Table 15.1 Comparison of predicted and computed value for wave- 
number spacing of adjacent modes. 
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We see that at = 5 . 8 , (15.19) is satisfied starting with the difference 
between the third and fourth modes; at = 10 it is satisfied approximately 
for the first three differences, and then exactly starting with the dif- 
ference between the fifth and sixth inodes. The upward -sloping line thus 
marks an important boundary. Above this line, the spacing solution is 
satisfied exactly; immediately adjacent to the line, it is not satisfied 
at all; and below the line it is satisfied approximately at high Mach 
numbers. 

With the eigenvalues of the multiple neutral solutions established, 
the next step is to examine the eigenfunctions. For this purpose the 
magnitude of the amplitude function 2 + C = tt / * W ) is sketched in 
Fig, 15*3 for the first 8 modes. The first thing to note from this figure 
is that the number of zeroes in is one less than the mode number u . 

For example, the second mode has one zero, and is l80° out of phase 

with ^ ; the third mode has two zeroes and is In phase with 

2^ . The number of zeroes in is the surest identification of 

the mode under consideration. Ejy keeping track of the phase relation 
between 2?^ M and 2^1^) , it is possible to determine when there is a 

change from one mode to another. 

Second, the appearance of the eigenfunctions in Fig. 15*3 tends to 
confirm what is suggested by the simple theory given above: there is an 

infinite sequence of wave lengths in the supersonic relative-flow region 
which can satisfy the boundary conditions. Third, the magnitude of z^Co'l 
is a minimum for the fourth mode. ( the same for all modes.) 

Since the fourth mode at - 10 is on the upward -sloping line of Fig. 15*2, 










this is another indication cf the special nature of the neutral solutions 
that make up this line. For other inodes, tends to become large 

away from* this line. In particular* this ratio becomes indefinitely large 
as n -*«) „ This same phenomenon will reappear in a slightly different 

form for three-dimensional disturbances and is illustrated in Fig. 15 . 13 . 

15*2.2 Singular neutral solutions 

A further consequence of a region of supersonic relative flow in the 
boundary layer is the existence of another class of neutral disturbances. 
Although these disturbances are also subsonic* the term neutral subsonic 
disturbances will always refer to the disturbances associated with the 
generalized inflection point. The new class of disturbances is Called 
singular because the solutions cannot be computed directly* but are 
obtained in the limit as C;~o from a family of amplified solutions. In 
the limit, c*. -*• l and ► co * That is* these neutral dis- 

turbances propagate with the free-stream velocity* and* since for a 
finite tr£o'i the free-stream pressure fluctuation Trl^A is zero* the 
disturbance is restricted to the boundary layer. The wave numbers of 
these solutions are designated by (X in . The first subscript refers to 
c-, =. 1 , the second is the mode number. The smallest ct, H is called the 
second mode* the next smallest the third mode* etc. 

Figure 15 • 4 gives i and av^irc-*^ as functions of 

V} for almost neutral disturbances of the second and third modes at 

~ 5*8. Both the amplitude ratio and a..e , n itud e a r the phase change 
are nearly independent of c*. near c h = 1* but the location of the phase 
change in the outer portion of the boundary layer varies someviiat as 


- 1 ^ 3 - 





















cv 1. This phase change is always located inside of the critical 
point. We see that for the second mode the pressure fluctuation under- 
goes only a 10° phase change, while for the third mode it undergoes a 
42° phase change in the outer portion of the boundary layer and a 180° 
phase change at y - 0,135, This 180° phase change is just what is 
encountered in the second -mode neutral subsonic solution, 

15 « 3 Amplified two-dimensional disturbances 

Nov that we are acquainted with the main classes of neutral dis- 
turbances, we can go on to the much more important amplified disturbances. 

i 

With the neutral solutions toiovn, it is easy to calculate amplified so- 
lutions by letting c ; increase from zero. It will turn out that not all 
of the neutral subsonic solutions lead to amplified solutions of interest, 
but at least one always does, as will be seen presently. 

It is convenient to plot the eigenvalues of the amplified disturbances 
in the form of diagrams of vs c- t and Ct. vs c* , Such diagrams are a 
necessity in carrying out the calculations. The complete picture of the 
inviscid solutions requires the damped solutions also, but to simplify 
the discussion, they will he deferred until Section 17*4, where the 
behavior of the viscous solutions is discussed. 

Figures 15 .5* 15 . 6 , and 15*7 give examples of eigenvalue diagrams at 
three Mach numbers, 3*8> 4.8, and 6,2, At = 3.8, the first two modes 
are included in the range of oc’s shown. If the calculations are started 
at both ot 3M c s and t c 4 , the two curves shown in Fig, 15«5 result. 

The curve which starts from (* 5l/ c s leads to the sonic neutral solution 
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a ~ o j c* ~ l-'/Mj, The solutions which lie on this curve have little 
phase shift in the pressure fluctuation and may he called first -mode 
amplified solutions* The curve which starts from d Si , Cj leads to 
the singular neutral solution o(i t , r . « l . The solutions which lie 
on this curve have a large phase shift in the pressure fluctuation except 
near the singular neutral solution. They may be called second -mode amplified 
solutions. 

It can be noted in these figures that for the first mode, c a <. c y < c 3 . 
This result is always found for the first mode. For the second mode, as 
seen from Pigs, 15*6 and 15«7j it i s possible to have <. c 3 . What is 
of major importance is that the amplification rate, <*C; , of the second 
mode is larger than for the first mode, a result which is always true. It 
is this feature that makes supersonic boundary-layer stability so 'different 
from incompressible flow. Pot only is there more than one mode of insta- 
bility, but it is one of the additional modes which is the most Unstable. 

At = 4.8, Pig. 15*6 shows that <* u is smaller than c< s , . This 
condition persists until = 6.0, after which a' marked change occurs as 
shown in Pig. 15*7 for = 6.2. The eigenvalue curve of the first mode 
joins the eigenvalue curve of the second mode at c ; ^ o . As a result, 
it is not possible to reach by following the eigenvalue curve .which 

starts at the neutral sonic solution, nor is it possible to reach the 
singular neutral solution by following the eigenvalue curve which starts 
at cC 5l . Instead, ot ^ is reached by starting at o( 5l . By reference to 
Fig. 15 * 2 , it can be seen that at = 6.2 oj Si lies on the upwaid- sloping 



portion of the o< 31 vs. curve. Although the first-mode neutral 
subsonic solution still exists, it has ceased to he of importance in 
the sense that there are no "adjacent’' amplified solution? of any 
consequence. At - 6.2, the magnitude of the ratio ir (a ) / for 

the neutral solution o( 5l f c s is l/8 of the magnitude of Tr^/Trc^ 
for the solution ot 6l ( c 3 .An obvious feature of Fig. 15*7 is that 
nowhere on the combined first- and second-mode curve does exceed c 3 „ 

As increases beyond 6.2, the. eigenvalue diagrams continue to 
change. At = 7.0, i.s less than o( sl , so the second- and third- 
mode eigenvalue curves cross. At ML - 7*5> the combined first- and 
second-mode curve does not lead to ot 3A l c 3 , but instead joins the 
third-mode eigenvalue curve at a c ( > o and leads to c< aa , c s . Hence, 
a single eigenvalue curve extends from the neutral sonic solution to the 
third -mode neutral subsonic solution. Neither of the neutral solutions 
ci S| t c s or , c 5 have "adjacent" amplified solutions of importance. 

It is also true that all along the combined curve, c„ * c 3 . 

When there are no longer separate eigenvalue curves for each mode, 
it can be difficult to identify a particular solution with a definite 
mode. However, the c<- c- diagram at =6.2 still retains the appearance 
of two separate modes. A close investigation shows that there is a rapid 
change in both the phase difference, ir ttCo 1 ! , and in the 

magnitude of the ratio tt-^ jvrcy^ near the "intersection" of two mode 
curves. These two quantities are shown in Fig. 15.8 as functions of c< 
at = 8, where the first three modes are merged. 

With the eigenvalues of the amplified solutions established, the temporal 
amplification rate etc. as a function of c* follows directly. When, at 
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several Mach numbers, oCc f is plotted as a function of oi for the first 

four modes, and the maximum amplification rate of each mode is read off, 

* / 

Fig. 15*9 can be constructed. This figure demonstrates that the second 

mode is always the most unstable mode. For > 6.5, the; first mode is 
not even the second most unstable mode, and at « 10, both the third 

and fourth modes are more unstable than the first mode. In order to get 
an idea of how these amplification rates compare to those found for in- 
compressible flow, we recall from Section 9*2 that for the Blasivs boundary 

-3 

layer at - 0, c t - 15 x 10 . The growth of a disturbance in a 

distance equal to one boundary-layer thickness is ' 



n 3 


( 15 . 20 ) 


where the group velocity has been replaced by the phase velocity. Table 

15*2 gives a few values of (ilA/A] 

M* s 


^ A / A L. 5 


0, 

finite R 

0.05 

0, 

£> = - 0.1988 

O.kO 

3, 

2nd mode 

0.03 

5, 

1st mode 

0.02 

5, 

2nd mode 

0,10 

10, 

2nd mode 

0,07 


Table 15*2 Amplitude growth over one bound ary-layer thickness. 


Thus for 3.5 <■ M^IO, a second~mode disturbance can grow as much over 
one boundary-layer thickness, according to the inviscia theory, as is 



Fig. 15*9 Effect of Mach number on maximum time rate of 
amplification of first four modes. 
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possible in the incompressible flat -plate boundary layer* It also 
follows from fig. 15 * 9 * that below about = 2.3 the boundary layer 
is for all practical purposes stable to two-dimensional inviscid dis- 
turbances. 

Of equal importance with the maximum amplification rate is the 
frequency of the disturbance which has this amplification rate. It 
is of no consequence to, have a large amplification rate for a particular 
frequency if the spectrum of an actual disturbance in the boundary layer 
does not include that frequency, figure 15.10 gives as a 

function of for the first four modes. The' 'dimensionless (circular) 
frequency is 


• ! i 

foj ' (15.21) 

lit R 


Thus for a fixed -frequency disturbance and a constant free stream, the 
dimensionless frequency varies as x Vl , and the disturbance can ex- 
perience the maximum amplification rate of Fig. 15*9 only at one par- 
ticular boundary-layer thickness. 

* . *, 1 

As a final observation from Fig. 15*9* we can note that the decrease 

1 : ■’ i 

in for M^> 4.5 is almost entirely a result of the increasing 

thickness of the boundary layer. If fc< s were plotted Instead of 

in Fig. 15*9* the amplification rate for each mode would be 
but slightly dependent upon past the peak values of Fig, 15 »9* 

Wow that we are familiar with the results of the inviscid theory for 
two-dimensional disturbances in the Insulated flat -plate boundary layer 
with wind-tunnel temperature conditions, we are in a position to understand 
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the changes that occur when any of these restrictions are lifted. Only 
flat-plate boundary layers have yet been investigated* but in the next 
three sections we shall find what happens when the free-stream temperature 
is different from that in the wind tunnel; when the disturbances are no 
longer two-dimensional, and when the wall is no longer iijpulated. Un- 

i 

fortunately, except for the numerical results there is a total absence 
of an inviscid theory for amplified disturbances. As general guides, it 
will be necessary to make use of two empirical Observations for the first 
mode, and for the higher modes, the result suggested by ( 15 . 19 )* that the 
wave numbers are inversely proportional to the* thickness of the supersonic 
relative-flow- region. ■ 

The two empirical results for the first mode are: first, is 

always between c„ and c s ; second, 'is related to the f difference 

c s - c a . That is, as c 6 — „ c Q ■ , (o( s c . ) o . We can see' from 
Fig. 15.1 that c 5 - c a has a minimum near 1 .6. This difference then 

increases and finally tends to level off. If for the 'first 

mode were to be plotted on a much larger scale than is used in Fig. 15 *9 j 
it would be found to have a similar behavior to*' c s - c 0 with a maximum 
near = 1 . 5 . This correspondence between c$ - c 0 and fa's c; 
suggests that anything that changes the relative positions of c a and c 0 
will also have an important effect on the amplification rate of 'the first 
mode. Both the free -stream temperature and the ratio of wall temperature 
to free -stream temperature affect c s . With three-dimensional' disturbances, 
c ti will change to c« 3 I - '/m, because the! effective Mach" number is 
the component in the direction of the wave nofmul. For a sufficiently 
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oblique wave c a is equal to zero. 

15*^ Effect of free -stream temperature on two-dimensional disturbances 

We shall study first the effect of changing the free -stream temperature 

while keeping the wall insulated. Figure 15.11 is a counterpart of Fig. 

u 

15.1. In this figure, there are three curves of c 3 (labelled A at > 7 * ) vs. 
for three values of T,* , ^0°K, l 6 o K and 300°K, Increasing "H* moves 

the c 6 curve closer to the c tt curve. Consequently, increasing T,* will 

« o 

stabilize the f st mode. Indeed, for "Ti = 300 K, c 0 between 

= 1.6 and 2 . 5 ^ and the first mode will be completely stabilized. 

In Fig. 15.12, the maximum amplification rates of the first and second 

T . <£. . 

, for three Mach numbers, 3 *5* 6.0 an( ^ 
8 . 0 . As expected from Fig, I 5 .ll, increasing. FT)* has a strong stabilizing 
influence on the first mode. The effect on the: second mode is. much smaller 
and can be either stabilizing or destabilizing,; depending on the Mach number. 
15 • 5 Three-dimensional disturbances 

15.5*1 .Neutral solutions ’ 

To treat three-dimensional disturbances ve, must use the tilde coor- 
dinate system of Section 13*3 to reduce the three-dimensional equations 
to two-dimensional equations. In this coordinate system, M, “ M, 

~ */co$ ^ and ^ is the angle between the wave normal and the free- 
stream direction. Since the boundary layers under consideration -are two- 
dimensional, the boundary-layer profiles are independent of \j/ . First 
we shall look at the variation of the neutral wave numbers, ol s ^ , with 
. Figure 15.13 gives the results at = 8 for the first three modes 
( \j/ is called <a* in this figure, and also in Fig. I 5 .I 9 ). We see that 
with increasing \jr , d S3 and c< Sl approach each other closely ; at ^ - 35 °> 
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1.0 



FREE-STREAM MACH NUMBER, M v 


Fig. 15.11 Mean boundary-layer velocity u at the gener- 
clised inflection point rj Sf and a!so 1 — U/Mi), 
vs free-sfream Mach number Mi 
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FREE-STREAM TEMPERATURE Tf 9 , °|< 

Wg. 15.12 Effect of free-sfream temperature T ^ on the 
maximum amplification rates of the first an,d 
second modes for M x = 3.5, 6.0, and 8.0 ' ' 
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DIMENSIONLESS NEUTRAL WAVE NUMBER 



Pig. i 5 ,i 3 ' Effect of weave angle on neutral , 
wave numbers at M x =* 8,0 
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and and 2 £l approach each other even more closely (to within 1$) 
at f = ^6.5°. , 

The slightly upward -sloping line that starts from at ^ => 0 
forms a sort of barrier for ■ When the o5 sl curve reaches the 
"barrier," it turns away from its previous direction to form a continu- 
ation of the upward -sloping line, while the « 33 curve forms a contin- 
uation of the ct £l curve. The same thing happens again when the ct £ , 
curve strikes against the "barrier" formed by the ot 52 upward -sloping 
line. At = 10, the first four modes form the same type of pattern. 
There, it is the curve starting from 2 54 at : '^■-c that is almost 
horizontal (more so than in Pig. 15*13) except for brief interruptions 
at y = 32 °, wnere it meets o( S3 ; at 47*2°, where it meets and 

at 60. 2°> where it meets oC 6 , . At V = 60.2^', of fil is within Ifegs than 
0*1$ of » ; v 

The various modes are easy to identify in. these figures by:. the number 
of 180° phase changes in the pressure fluctuation. Even though the eigen- 
values almost coincide at the angles listed above, the solutions do not. 

In some limit, probably c 3 -**i as , the; eigenvalues will be truly 

degenerate, i.e. , two distinct solutions will correspond to ohe:: eigenvalue. 
The wave-number curves will actually intersect, but the mode identity of 
a given curve will still change beyond the intersection point,. On 

physical grounds, it is necessary for the first-mode wave number to 

o 

always attain a finite value as y 90 > and for all higher-mode wave 
numbers to go to infinity, with a common asymptote at the angle where 
the supersonic relative-flow region disappears. :i 
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It is of interest to note that Fig. 15*13 hears a resemblance to 
Fig. 15. 2v There the upward -sloping line' plays the same role as the 
more nearly horizontal line of Fig. 15*13* If the angles of closest 
approach of Fig. 15*13 converted into component Mach numbers, we 

can construct Table 15-3* It is apparent from the table that the 



Mj_ = 5-8 

8.0 

•10.0 

2-D 

<*• ** 

^ = 3.9 

4.4 

5*0 

\ - 4.7 



6.5 

6.8 

6.7 

* 5 * * 



8*5 

,i >: 

8.6 


• Table 15*3 Mach numbers in direction of wave normal at whiqh 

wave numb' _ of adjacent modes are most nearly equal . 

■ . , ;; •> t 

phenomenon is the same fox Doth two- and three-dimensional disturbances. 
Whatever the mechanism is that produces these results, it must depend 
primarily on the Mach number and very little oh the bound ary -layer 
profiles. 

1 

The same two features that were noted for the upward -sloping line 
of Fig. 15.2 are also found for the related line of Fig. I5.I3, First, 
it is the mode which is on this line for which the ratio rc [o] f rr 
has its minimum value for any particular wave angle, as' is shown in 
Fig. 15*14. Second, it is the wave number of this line which is the 
end point of the eigenvalue curve of amplified • solutions which 1 starts 
at the neutral sonic solution. This behavior, which for two-dimensional 
disturbances establishes itself clearly only at high Mach numbers, is 


Fig. 15. 14 Amplitude ratio of surface pressure fluctuation a 
function of wave angle for first three, 'modes at = 8. 
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veil defined for three-dimensional disturbances as vill be seen in 
the next section, 

15 • 5 * 2 Amplified solutions 

With the neutral subsonic solutions established, the next step is 
to consider the amplified solutions. As y increases from 0 *, c 0 « l - 7 m, 
decreases, but since the mean velocity and temperature profiles remain 
fixed, so will c* . Consequently, the difference c s c t increases 
and ve can expect a destabilisation of the first mode. At the same time, 
the thickness of the supersonic relative-flow region will decrease with 

V 

9 and we shall not he surprised to find a stabilization of the second 
and higher modes. 

Figure 15.15 shows the time rate of amplification, 0 <c f f of the 
first and second modes at — 4. 5 as a function of the dimensionless 
frequency, c h * for several wave angles. The three-dimensional first- 
mode disturbances are indeed more unstable than the two-dimensional 
disturbances, and the second -mode three -dimens ional disturbances are 
more stable than the corresponding two-dimensional disturbances. Con- 
sequently, the most unstable second-mode disturbance is always two- 

dimensional. This same result holds for all of the higher modes. The 

o 

most unstable first-mode disturbance is at an angle of close to 60 , 
with an amplification rate twice the maximum two-dimensional rate and 
with a frequency a little over one -half of the frequency of the most 
unstable two-dimensional disturbance. The maximum in at a value of 

, o o 

Y between 0 and 90 comes about from the following circumstances. As 
C s - c 0 increases, the maximum value of c f also increases. At first this 









> 

increase leads to a corresponding increase in . However* since 

<* (not Si ) must go to zero as if" 'l 0 ** the decrease in oc eventually 

outweighs the increase in c f and the amplification rate starts to decrease 

with further increases of y and is zero at Y = $ 0 ° . 

At = 4.5* the first two modes are completely separate at all wave 

angles. At a 8 * the first three modes are merged for two-dimensional 

disturbances. Figure 15.16 shows what happens as Y increases. At 

Y - 30°* the first three modes are still merged. However* at Y - 45°* 

only the first two modes are merged. This result is in accord with 

Fig. 15.13* where the second -mode wave number* 0i Sl * is on the upward- 

sloping line at Y = 45° and thus is the end point of the eigenvalue curve 

, 0 

for the merged amplified solutions. At Y = 5^ * the situation is 1 the same. 

The next angle* Y = 56 , is of interest because this angle is very near 

«-<■ 

the angle of closest approach for o( 4( end d 3t . The end point should be 

o( sl according to Fig. 15.13* and Fig. 15.16 confirms that it is. '■ However* 

the tfCj vs. e* curve in Fig. 15.16 gives no obvious evidence of the second 

mode as the local maximum in otc ( that has previously always marked the 

second mode is no longer in evidence. Hence we have an example Where it 

is necessary to examine the phase of the pressure, fluctuations ih order to 

make a positive identification of the mode. Figure 15-17 gives the ratio 

TT CO * aricl Fig. I 5 .I 8 the phase difference ai-^-rrCo} - * 

as functions of oc „ These figures show that there is a rapid change from 

a first-mode disturbance to a second-mode disturbance near - 0 ;l 6 even 

o 

thougii acj continues to decline monotonically. At the next angle* Y = 60 * 
Fig. 15 . 13 indicates that ot 5 , will be the end point of the amplified region 
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■which will consist entirely of first-mode disturbances. Figure 15. l6 
and the phase of ‘HTo} confirm this to be true. The angle at which 
the supersonic relative -flow region disappears, and with it the higher 
modes, is just over 60°. • . 

Once we have found oUj as a function of etc,. for several values 
of y at a particular Mach number, we can construct Fig. 15.19* In 
this figure, the maximum temporal amplification rate of the first two 
modes is plotted against y for the four Mach numbers 4-, 5# 5*8, 8.0 and 
10.0. At all of these Mach numbers the most unstable first-mode dis- 
turbance is at an angle of between 50° ancl 60°/ and has a maximum 'ampli- 
fication rate .roughly double the amplification rate of the most Unstable 
two-dimensional disturbance. The two-dimensional second-mode disturbance 
is the most unstable at all of these Mach numbers, and the decrease of 
with increasing y is sharpest at thenlower Mach numb&rfe. 

A summary plot of the first-mode maximum amplification rates appears 
in Fig. 15*20. The wave angle of the most unstable disturbance (to within 
5°) is shown in the figure, and the maximum tvo*ddmensionol amplification 
rates are shown for comparison. An interesting >. change in the relationship 
between the 1 two- and three-dimensional amplification rates takes place for 
M- l < h. Above this Mach number, (ctCjW* of the i three -dimensional dis- 

l 

turbances is about double the two-dimensional amplification rate.! sun 
contrast, at -M^ = 3«0 the ratio of the two amplification rates is 5*»8; 

« 2,2 it is 33; and at = 1,8 it is 130. ' We recall from 'Fijg. 15.I 

that it is near M^ » 1,6 that the difference c a >- c 0 is the smallest. 
Therefore, the sonic limit acts as a severe constraint on the two-dimensional 
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Fig. 15-20 Effect of Mach number on maximum amplification rates of two- and three- 
dimensional first-mode disturbances, Inviscid theory j insulated wall. 
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disturbances in the low Mach number range. Indeed the two-dimensional 
disturbances are almost completely stable even though the generalised 
inflection point is well out in the boundary layer. When this constraint 
is removed 3 as it is for three-dimensional disturbances, the amplification 
rates increase sharply. We may consider the three-dimensional maximum 
amplification rate as the one that most clearly reflects the inherent 
instability of a given boundary-layer profile. 

The instability of a compressible flat -plate boundary layer! which is 
due to the generalized inflection point is a weaker instability than the 
instability of an incompressible boundary-layer profile with an Inflection 
point. In Table 15.4 we compare the growth ove'r one boundary -layer thickness 
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r ■ 

o, § - - o.l 4 

* ,f7 0.10 

* { 
1. 

o, ^ - 0.1988 

■ 0.40 

n c 

1.6 

0 

0 

9 

0 

v r 

2.2 

0.012 


3.0 

’ 0.022 

i j 

4.5 

- >: 0.036 

*. n 

10.0 

1 : 0.042 

T'l 


Table 15,4 Comparison of growth of most unstable first-mode irivlscid 
disturbances in insulated flat -plaice boundary layers and 
in Falkner-Skan boundary layers. 

. ! ■ 'I 

of the most unstable three-dimensional first-mode disturbances at a few 
Mach numbers with the growth of the most unstable inviscid disturbances 

> i i 

in two of the Falkner-Skan boundary layers* We see that first -mode dis- 
turbances in the compressible boundary layer grow much slower than in the 
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two Falkner-Skan boundary layers. Indeed the maximum growth rate is 
smaller than the maximum growth rate of 0.05 (at finite Reynolds 

i 

numbers) in the incompressible Blasius boundary layer. 

15.6 Effect of wall cooling 

4 ' 1 ' i , 

15.£.1 Two-dimensional disturbances 1 

, h 

Perhaps the most celebrated result of the stability theory for 
compressible boundary layers was the prediction by Lees (l$?Vf) ihat 
cooling the wall stabilizes the boundary layer* This prediction was 
made on the basis of the asymptotic theory , dnd a criterion was provided 
whereby the ratio of wall temperature to recovery temperature ?at which 

' . j 

the critical Reynolds number becomes infinite 'can be computed J* 1 Although 
Lees's original calculations contained numerical errors, the temperature 
ratio for complete stabilization was later computed correctly by a great 
many authors. The most accurate of these calculations gave the result 

; . ' 1 ' f J • » 

that complete stability can be achieved for rlr< M^< $ by sufficient 
cooling. I These calculations can he criticized in two important respects. 
First, no indication is given as to how the raraplif ication rate Ovaries 
with wall temperature. Second, and most impbrtant, no account is taken 
of the existence of the higher modes. In this section, we shall see that 
the invlscid theory can remedy both of these deficiencies. ; 

> The first step is to find out how the mean boundary layers changes 

with wall cooling. In Fig. 15*21, the stability function v 0 which 

} . : 

is related to (li'/T) by ( 15 • 5 is giveri nas a function y at 

1 ■ ■ 

•M^ = 5*8 for several values of 0 W , where is the enthalpy difference 
ratio -Ail / Wo - A,) * It is, < given in terms of $ome common 

1 ' ^-GfNAL page 13 .. 
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INVISC1D STABILITY FUNCTION, 



i’ig, 15=21 Distribution fiforovgh boundary layes* of 
stability function v 0 for foyr values of $^\ 


*Si SXi# 
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temperature ratios at = 5.8 in Table 15.5*' In this table ^ is 

lh. ’* 1- 


Bw 

T*/T„ 

■ ; Tv,/T, 

‘ 0.847 

1.00 

■n. 6.7 

0.50 

0.65 

"1 4.4 

0.10 

0.25 

!.!!• 1*7 

0.05. 

0.20 

‘ 1,3 

0.01 

0.16 

1.1 

0 

0.15 

i » 1.0 

-0.10 

0.05 

-■! 0.32 


Hi- 

Table 15.5 Temperature ratios as functions of & w at *p $.Q 


the recovery temperature, i.e., the equilibria or adiabatic vail 
temperature. In Fig. 15*21 ve see that for the insulated vail, f Z^'/T ^ 
has a single aero located at . For 0.01* Sw < 0.847 j it has tvo 
zeroes located at and > >> J2 )• For 0 W < 0.01, 

Oi'/T V has no zeroes. The zero at is alvays belov , the point 
vhere u s c 0 » l - 1 /n\i • The zero at '• is above ^ 0 provided 


> 0.05. 


1 1 I., i 1 ■ 1 . 1 


Next ve look at the relation of c 5 , considered to be the mean 


velocity at , to c 0 . As long as > there can be no 

neutral disturbance connected vith « It .-cannot be subsonib because 

c 5l < 1- Vim, . it cannot be supersonic because, by (15*4), the 


Reynolds' stress would be zero throughout the boundary layer, a Condition 

( 

incompatible vith a neutral supersonic disturbance as ve shall ’bee in 
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Section 15.6,3* Both c a and c a a re shown in Fig. 15.22, where c 5 is 
plotted as a function of T W /7P at « 5.8, 8.0 and 10.0 for 

T t = 50 K, and at «= 8.0 for T t * = 80°K. We see that at - 5 » 8 > 

c s is reduced "below c fl with sufficient cooling, and the first mode can 
"be completely stabilized. However, at - 8 and = l<p, c 5 remains 
greater than c 0 and the first mode cannot be completely stabilized no 
matter how much the wall is cooled. 

In order to -find out what happens to the second mode as the wall is 
cooled, we must look at some eigenvalue diagrams. In Fig. 15*20, * is 
plotted against Cj for five values of d w at S = 5.8. In Fig. 15*24 
three of the corresponding c r vs. c : diagrams 'are given with the two 
remaining diagrams in Fig. 15.25. Only amplified and neutral solutions 
are shown. From these diagrams we see that the ( first mode disappears 
when y 5l moves below y 0 and there can no longer* be neutral subsonic so- 
lutions, ■ It seems to be a requirement that ar family of amplified solutions 
can existf-only if their eigenvalue curve starts at one neutral solution and 
ends at another. When *: y 0 , or does f not exist, there 1 is still 
an eigenvalue curve which starts at the neutral sonic solution,.; but the 
solutions; along this curve are all damped as we shall see in Se&bion 17. 
When somd of the modes are merged for the insulated -wall case,, as happens 
for > 6.0, cooling will separate the modes' so that there are Tdis tine t 
eigenvalue curves for each mode just as at lower Mach numbers rfqr the 
insulated-wall case. -i : R n 

The second mode always has the singular 'neutral solution a si one of 
its necessary two neutral solutions, just as ^fcfoe first mode always has 

■! - -? • r, h <: : 
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Fig. 15.22 Mean boundary-layer velocity at generalized 
inflection point as function of ratio of wall 
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the neutral sonic solution. These solutions exist regardless of the 

temperature boundary condition. All that is required for the singular 

i ; :i ' ‘ I 

neutral solution is a region of supersonic relative flow fop - 1. 

As the wall is cooled and c s becomes smaller., the difference 1 - c 5 
increases and there is no tendency to reduce thO maximum value of !c; 

l . j 

• for the second mode as there is for the first mode as a result of the 

! 

• I 

decrease in c a - c 0 . Indeed we see from Fig. 15*24 that initially C; 

I ‘ ! 

increases. Since «, also increases, there is an increase of the maximum 
amplification rate of the second mode with wall cooling. The increase 
in the wave number at tecjW** is almost entire!^ a result of the- 1 ' change 
in the boundary-layer thickness with cooling. ]>; f.\ 

, I 

When moves below >j a and' then disappears, it is no longelr 

possible for the second. mode eigenvalue curve toxhave as an end 

point, since the neutral subsonic solutions no longer exist. Nor is 
there any possibility of the neutral sonic solution providing the tnec- 
essary end point as this solution has only firdti-mode solutions )i $ its 

neighborhood. If we look at Figs. 15-23 and 15-24, we see that for 

V 

si <M0 r: , where ^ 5t > ^ , there has appeared a new neutral 

solution together with a family of amplified solutions.' For all'' Of 

• „ ■ i i ■ ■■ i j • • i ‘ . 

* thdse solutions, is less the'i l - Vm» » and; is supersonic with : respect 
to the freq stream. Consequently, we have a neutral supersonic (Solution 

| \ . * 4 I 1 1 W I L i | * I 

and supersonic amplified solutions of the type^dlscussed in Section 13-7-2. 


All of these solutions are outgoing waves. It ; '-£s evident from Figs. 15*23 

(I H .... I . I 

and 15.25 that when <. y o or does not exfi^t, the eigenvalue? : curve 

>1 . . i . . iii. >1 ii. 

of the amplified subsonic solutions joins the Eigenvalue curse of* 1 the 

r , 1 1 it 1 1 , 1 1 1 . i 1 / . 
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1 > j 11 . 


• I , ' ! I I I • ■ l l • • . . * ( i i 

.1 ' ■ i i - « -i i ! i ' I ■ l i ■ , i i i i , i. 

ill i 1 1 i . , • ■ ( i ■■ i I . i ■ . . . 

amplified supersonic solutions at a c ( - > o . The end point of this 

’ i i . i * i . i 14- ■ 

combined eigenvalue curve is another neutral supersonic solution which 
has not yet been studied, but it appears to bbof a singular nature. 

•y J, 

The important result is that conditions exist by which tjiere can be 

i l > 4 

second -mode amplified solutions regardless of .whether or po ( t there is 
a generalized inflection point in the boundary layer. Cobiing the wall 

. • / : ‘i 

does nothing to remove this source of amplification. ’ This same result, 

. . i • ' • / 

of course, applies to all of the higher modes. ‘ ' 1 • 1 11 l ' t ' ** 

If the maximum amplification rate of the first three modes is 
computed for several 77, / "TV at = $.8, then Fig. 15.26 can be con- 
structed. We see from this figure that as 71, / Tv- decreases, 
of the first mode decreases to zero monotonically. In contrast, 

of ^he second mode increases by 50$> f and of the 

third mode increases by a factor of three. The amplification rate otc ( - 
is the proper one when we wish to make comparisons for a fixed free stream 
and at a fixed R K . If we wish to compare disturbance growths over a 
boundary-layer thickness, then we need c- . Table 15.6 gives some 


0W 

Tw/TV 

% 

% 


0.847 

1.00 

20.0 

9^ 

8.5 

0.50 

0.65 

16.1 


8.3 

0.25 

0.40 

13.2 

6.0 

7.7 

0.10 

0.25 

11.2 


• 7.0 

- 0.10 

0.05 

8.5 

3.2 

5A 


Table 1J5.6 Maximum second -mode amplification rate based on boundary- 

layer thickness for cooled -wall boundary layers at ML .» 5*6* 
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MAXIMUM TIME RATE OF AMPLIFICATION, X 10 



Fig* 15.26 Maximum time rate of amplification of first 
modes at M 1 = SM as function of rati© 
of wall to recovery temperature 
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values of C; which correspond to the values of C«*C;W, fcA 

of Fig. 15*26. Also included in the table are some values of , 
the thickness of the supersonic relative-flow region for the phase 
velocity of the most unstable disturbance. We see that decreases 
a bit more rapidly than ^ as T w / TV decreases. 

As two final results concerning the effect of cooling on two- 
dimensional disturbances,, we may look at Figs. 15*27 and 15. 28, where 
<*C; is plotted against c* at = 8 and ~ 10 for the insulated 
wall and for T w /Ti- =0.05* At = 8.0 and = $0 K, the first 
three modes are merged for the insulated wall; at = 10 the first 
four modes are merged* In the highly cooled case, the modes are separate 
at both Mach numbers. The first mode is slightly unstable at = 8.0, 
but the amplification rates are too small to show in the figure* At 

= 10.0, the most unstable first-mode disturbance has an amplification 
rate equal to l/3 of its value for the insulated -wall case and a wave 
number which is larger for the cooled wall than for the insulated wall. 

At = 5*8; this wave number decreases with increased cooling. The 
higher modes at both Mach numbers are destabilized and shifted to higher 
wave numbers, just as at ^ 5*8. The maximum amplification rate of 
each higher mode is approximately doubled as a result of the wall cooling. 
Since the boundary-layer thickness is also cut in half, (0^ c,\ Aney> is 
virtually unchanged by the cooling. 

15.6,2 Three-dimensional disturbances 

Cooling the wall stabilizes the first-mode two-dimensional disturbances 
by moving vj Si toward y o . Since a three-dimensional disturbance at a 
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Fig. 15.27 Time rate of amplification vs wave number for 
first three modes for insulated wall and 

e 



TIME RATE OF AMPLIFICATION, oc/XIO 3 




sufficiently large wave angle ^ reduces ^ well below and makes 
it impossible for y St to be less than , we may ask if 'in this case 
there also exists a stabilizing mechanism for first -mode three-dimensional 

disturbances. The answer is yes, and it involves the second generalized 

f 

inflection point, y JJr , which we found to have no significance for two- 
dimensional disturbances since it was always less uuan y 0 * However, 
with increasing y , y 0 can be made low enough 'so that y st > y o „ 

As soon as this condition is fulfilled, there will be two neutral sub- 
sonic solutions, and the eigenvalue curve of the first-mode amplified 
solutions will be located between these two solutions . As T*, / T k 
decreases, y SL and the maximum amplification rate will decrease, 

just as it does for two-dimensional disturbances when y St —*■ y o « At 
as we see from Fig. 15.21, with sufficient cooling • 

Consequently, the first mode is completely stabilized even for three- 
dimensional disturbances, but at a lover value of Tv /“TV than was 
required to completely stabilize two-dimensional disturbances. At higher 
Mach numbers, y si is greater than regardless of the amount of cooling, 

and there will always be amplified three-dimensional disturbances. However, 
the amplification rates of these disturbances will be much less than the 
amplification rates of the insulated-wall boundary layer. 

For second-mode three-dimensional disturbances, the same result holds 
as for the insulated -wall case: the most unstable disturbance is two dim- 
ensional. Complete calculations have not yet been carried out for three- 
dimensional disturbances in cooied-wall boundary layers, but there has been 
enough done to establish the above described behavior and also to provide a 
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curve of vs. Tv/TV at ^ - 5*8. (The maximum here is 

with respect to both <* and .) This curve is included in Fig, 19.4. 
15.6.3 Neutral supersonic solutions 

As already mentioned, a neutral supersonic solution Appears in Figs. 

15.23 and 15.24 for £ w = 0.10 { T w /Ti- = 0.25). This solution fits the 

description given in Section 13.7*2. It is a purely outgoing wave, and 

in the coordinate system fixed relative to the -phase velocity it is a 

Mach wave. Since the wave extends to infinity with undiminished amplitude, 

it transports energy from the. boundary layer to infinity. However, it is 

a neutral disturbance and the net energy transfer must be zero. The energy 

carried to infinity is made up by energy production in the boundary layer 

due to a Reynolds stress which exists in the region ij c <. ^ 

This Reynolds stress, as calculated by numerical integration, is shown 

In Fig. 15.2Q. This figure agrees with the sketch given by lees and Lin 

(1946) for a neutral supersonic disturbance. As noted by Lees and Lin, 

must be less than c 5 in order that the Reynolds stress, which is 

given by (15-4), is positive. A neutral supersonic solution of this 

type has also been found at Kj = 5.8 for f? lv = 0.15, but not for : ~ O. 25 . 

A plot of the magnitude of the pressure amplitude function, ]'2^| , 

is given in Fig. 15-30 for the neutral supersonic solution of Fig. 15-29 

and also for the second -mode neutral subsonic solution. The two amplitude 

n 0 

functions are similar in character, and, because of the single 180 phase 
shift, we can assign the neutral supersonic solution to the second mode. 
This result is a reasonable one, in view of the: fact that the second -mode 

4 % 

or has disappeared. 


amplified solutions join the supersonic amplified solutions when> 
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We may now ask if there are any neutral supersonic solutions which 
are purely incoming waves instead of the outgoing waves just discussed. 

Wo solutions of this type have yet been calculated, hut inspection of 
complete eigenvalue diagrams which include the damped solutions makes 
their existence probable. Diagrams of this sort for Tw/T h = 0.25 
'and 0,0p ( = 0.10 and - 0.10) at = 5.8 are included in the figures 

of Section 19*2. Prom these diagrams we see that when the eigenvalue curve 
of the second -mode amplified solutions for - 0.10 is extended into 

the damped region, it joins up with the eigenvalue curve of a family of 
damped supersonic solutions. According to the results obtained in Section 
13.7*2, these solutions are incoming waves. At this value of & ^ , there 

is a neutral supersonic solution and a family of amplified supersonic 
solutions, all of which are outgoing waves. At 6 ^ - - 0.10, the eigen- 
value curve of the second-mode amplified solutions connects up with the 
previously separate amplified supersonic solutions. By analogy' with the 
situation at 8 *, ~ 0.10, we may speculate that a neutral supersonic 
solution also exists here, but that it is a purely incoming waveband is 
associated with a family of damped incoming waves. Obviously, this neutral 
disturbance receives energy from infinity, and the energy gain must be 
balanced by an energy loss due to a negative Reynolds stress in the region 
7$ 7 * 75 .In order to have a negative Reynolds stress, 

'according to (15.^), c*- must be greater than c 5 „ 
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16. Review of Asymptotic Theories J ' 

The inviscid theory as presented in the preceding section has J 
given us a large body of numerical result s, the 'applicability of which 

i 

has not yet been established. To find out where these results can be 
used, it is necessary to turn to the viscous theory. Historically, 

J 

the viscous theory was first attacked by the same methods that had 
proved successful for incompressible flow. A brief account of this 

theory and its Im mediate successors is given in this section. Only 

v ■■ ■ 1 

two-dimensional disturbances are considered, 

16.1 Lees -Lin theory 

Each form of the asymptotic theory depends upon reducing the 
complete stability equations of Section 13*5 to a simplified system 
of equations by means of order-of -magnitude arguments. The inviscid 
equation provides two of the necessary six solutions. The simplified 
equations are derived to be valid in a viscous region, and they will 
provide the four additional solutions needed to. solve the eigenvalue 
problem. The region of validity of the Lees-Lin equations is the 
Immediate neighborhood of the czitical point. 

A new independent variable $ is defined, where 

5 = 

£ 

and 

, 3 

£ - R \ 


(16.1) 


( 16 . 2 ) 
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This variable is essentially the same as the one used in the incom 


pressible theory in Section 6.4. 
lead to 

z, «' x, 

= X* 



The order-of -magnitude considerations 

(16.3; 

(16.4) 

(16.5) 

(16.6) 

(16.7) 


£2 S = 


(16.8) 


where the 'X^ are all of the same order and may be considered to be the 
first terms in a power-series expansion in £ . When these relations 
are substituted into the complete equations* the rean-flov quantities 
hr c and T expanded in series about y c * and only the lowest order 
terms kept* the equations of Lees and Lin are obtained. They are 

X. “ Xa. (16.9) 

(i s *Xi Xj) ±z*~ %*. (16.10) 

isi js, 

=* - i Xi (16.11) 
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0 


(16.12) 


X - 


X - X, 


( 16 . 13 ) 


X , ^(iu'.sXf * rjy 3 ) 


( 16 . 14 ) 


If (l 6 . 10 ) is differentiated, the following equation is obtained 
for X3 


X' _ mi i x 


=3 0 


(16.15) 


Since this equation is the same as (6.29), two viscous solutions are 
provided by (6.32) and ( 6 . 33 )- From (l 6 .l 4 ), the equation for X$* is 




t 


K 




XU X: 


(16.16) 


Consequently, X3 is coupled to Xy , and velocity fluctuations will 
bring about temperature fluctuations. The important point of the Lees- 
Lin theory is that temperature fluctuations do not similarly bring about 
velocity fluctuations ( Xy does not occur in the equations for X, and 
). The velocity boundary conditions at the wall may be satisfied 
by the inviscid solution and the proper viscous solution of (16.15). 

Then solutions of the homogeneous part of (l 6 . 1 o) are simply added to 
the temperature fluctuations associated with these solutions to satisfy 
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the temperature boundary condition at the vail. The temperature 
fluctuations do not enter the eigenvalue problem, which is identical 
in all essential respects to the eigenvalue problem in inqompressible 
flow. 


This form of the theory was used by Lees (j. 9^7) to compute neutral 
stability curves at = 0.5,, 0.7, 0.9, 1-1 and 1.3 for an insulated 
wall, and at = 0.7 for three cooled -wall cases and one heated*-wall 
case. In a later report, Lees (1952), he computed amplification rates 
and the overall growth of constant-frequency disturbances at = 0.7 
for the insulated-wall case. 

The Lees-Lin equations can also be written without expanding' 1i 
and T in power series. In this fo:rm, the equations for T , and & 
are as follows: 


f 


Hi 


v 


Q 


(16.17) 


i 


-i-f 


(16.18) 


0 


11 


is 


chRstT 1 

if 




( 16 . 19 ) 


where the independent variable is yj , The validity of these equations 
still rests on the supposition that U-c is a small quantity, however, 
we may also do as in incompressible flow and simply regard these equation 
as an approximation that is useful when If - e is not small, which is 
the way the theory is applied in any case. A further step with 'll - c. 
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not considered small is to obtain the counterparts of the Tol'JLmien 
viscous solutions of Section 6 . 5 . In this form, the theory is an 
approximation to the Dunn-Lin theory which is discussed ifi the next 
section. Neutral- stability curves have been computed by Mack (i 960 ) 

at = 1.3, 1.6 and 2.2 from this version of the theory. The inviscid 

* 

solutions were calculated hy the same numerical method as used by Lees 
( 19 ^ 7 ) (power-series expansion in oi* ), except that one additional term 
was used, 

16.2 Dunn-Lin theory 

Because the; viscous solutions of the Lees-Lin theory are strictly 
valid only near the critical point, Dunn and Lin (1955) developed an 
improved theory where the viscous solutions are of the Tollmlen type 
and are valid everywhere . The viscous equations were derived by 'Dunn 
( 1953 ) from careful order-of -magnitude arguments applicable in the wall 
viscous region. In this viscous region, the ordering parameter ;ts 
£ =5 f the same quantity which appeared in Section 6,2. 

With li- c no longer considered a small quantity, the viscous equations 
are 


f m i ctEUZf.-O ^ 


1 / 


(16.20) 


q' = - 1 i 


I (U - c) 

T 


(16.21) 


0 )l I et R V t cl — 0 

V 


(16.22) 
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When U-c is small , (lo.2l) reduces to (l6.l8), and when h- c is not 
small, ( 16 . 19 ) reduces to (lo.22). Equation (l6.20) is the same in 
both the Lees -Lin and Dunn-Lin theories. 

There is an important difference between the Dunn-Lin and the Lees- 
Lin viscous equations. In the Lees-Lin theory, velocity fluctuations 
induce temperature fluctuations, but temperature fluctuations do not 
induce velocity fluctuations. In the Dunn-Lin theory the contrary is 
true. Velocity fluctuations do not induce temperature fluctuations, but 
temperature fluctuations induce velocity fluctuations. In other words, 
the Lees-Lin equations have the four viscous solutions 


) $3 y ) / (^4 / ^4 1 

(0 , 0 , S f ) , (0,0, 6 i\ 

end the Dunn-Lin equations have the solutions 

(L - $>, 0 } , ft* , ?* / 0 ) 
(0 j > (0 , <j t , eh 


(16.23) 


(l6. 24) 


Consequently, in the Dunn-Lin theory the velocity and temperature boundary 
conditions must be satisfied simultaneously, and as a result the temperature 
fluctuations enter into the eigenvalue problem. 
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If we designate the inviscid solution that satisfies the boundary 

conditions at infinity by ( F ,$,©), and eliminate the viscous 

f 

solutions which become large as y *-» °o (subscripts four ^nd six 
according to the conventional notation), we can write the' boundary 
conditions at = o as 


Fw +- 

A 'fgw ■+■ 

B fjv = 0 


5w * 

^ ^3w * 

B ^ = 0 

( 16 . 25 ) 

©V + 

A 0 3w -* 

0 Bf* - 0 



In order to have a non -trivial solution of these homogeneous equations 
for A and B, the determinant must be zero. With and $$ zero from 
the Dunn-Lin theory, the expansion of the determinant gives 


Jw 


+»w few Fw 


(l 6 . 2 . 6 ) 


The inviscid equations given in Section 13.7*1 can be manipulated to 
yield 


6 w 

Fw 


{ y- i\ c Mj 4 


Jr-0 M? Ul - 



i Sw 

Fw 


(16,27) 


With (16,27) substituted into (16.26), the eigenvalue equation is found 
to be 


tl?. 

* jw d $vj 


t - 


CY~t] /Vi* Idw “ 




( 16 , 28 ) 
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This equation was first obtained by Re shot ko (i960), but another 
equation which reduces to the above for ~Tj - 0 vas derived earlier 
by Dunn and Lin (private communication). The inviscid temp are on 
the left-hand side^ and the viscous terms are on the right-hand side 
just as for incompressible flow. In the Lees-Lin theory, (16.28) is 
replaced by 


SL __ 

F'-v -f 3w 


O6.29) 


which is identical to (6.2) in the incompressible theory. 

When the temperature fluctuations are unimportant in the sense 
that the associated velocity fluctuations given by (p^- are small, 

(l6.28) also reduces to (16.29). For the insulated -wall case, the 
neutral-stability curve is not influenced by the temperature fluctuations 
up to about = 1.6. At - 2.2, the neutral- stability curve with 
temperature fluctuations included is outside of the neutral curve 
obtained from (16.29). Unfortunately, the neutral -stability curve 
obtained from numerical integration of the complete stability equations 
is inside of the neutral curve obtained from (16.29) as will be shown 
in Section 17.1. 

16.3 Lees-Reshotko theory 

A further attempt to improve the asymptotic theory was made by Lees 
and Reshotko (1962), reported first by Reshotko (i960). The separation 
into inviscid and viscous solutions vas retained, but all solutions were 
calculated by numerical integration of the appropriate equations. Because 
the viscous solutions were to be obtained by numerical integration, a 
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more complete system of viscous equations could be used than in the 
normal asymptotic method. The viscous equations, derived from order- 
of -magnitude arguments, are 


i m + L iff' f" » — it- it s“ 

/* IT /* AT 


U R ( lc~ O /V* 4 _ 4 ' _ ' e\ «- o 

* £ It t ; 


(16.30) 


+• H f - I tf . i ( 1j- c) ^ , 

T ^ T 


(16.31) 


it 

/* *IT 


0 1 ' + h tt. T' 0* + m, 26* -f 


i ct. R v ( ^ _ A. R T T gj 


(16.32) 


A neutral-stability curve was computed by Lees and Reshotko at 
= 2.2. It agreed vith the Lunn-Lin theory on the upper branch, 
and on the lover branch vas shifted somewhat to the right. Hence 
there was an improvement over the Dunn-Lin theory, although a small 
one. An attempt vas also made to compute neutral -stability curves 
at = 3-2 and 5 *6, but at these Mach numbers the neutral curve vas 
found to break up into multiple loops of a puzzling nature. These 
loops have not been explained, but it does appear that they are funda- 
mentally different from the multiple loops to be discussed in the next 
section, and which can be identified vith the multiple Inviscid modes 
already discussed. 
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17* Viscous Theory - Two-Dimen s ional Disturbances, Insulated Wall 

We now take up the viscous theory to the extent it can be developed 
from the results obtained from numerical integration of the complete 
stability equations* With these results, we can assess the accuracy 
of the asymptotic theories ana also find under which circumstances it 
is possible to make use of the inviscid theory. Further,, we can study 
such things as the influence of the Reynolds number on amplification 
rates and frequencies of the most unstable disturbances. Finally, we 
can compute the growth of disturbances of constant frequency as they 
travel through the boundary layer, and thus be in a position to judge 
the response of the boundary layer to any particular disturbance spectrum. 
17*1 Neutral- stability curves 
17,1.1 Comparison with asymptotic theory 

Figure 17*1 gives a comparison of three neutral-stability curves 
for two-dimensional disturbances at ~ 2.2. The dimensionless frequency, 
Co* 3/* / U* s is plotted against R = ft* ( Cu* is referred to as ^ 
in the figure). The outermost curve is from Mack (i960) and was obtained 
from (16.28), the eigenvalue equation of the Dunn -Lin theory which includes 
the temperature fluctuations. The innermost curve is from the numerical 
integration of the equations of Section 13*5 according to the method of 
Mack (1965a). There is seen to be a large discrepancy in the two results. 
Somewhat better agreement is found by omitting the temperature fluctuations. 
The intermediate curve was obtained by direct numerical integration of the 
simplified system of equations given by Dunn and Lin (1955)* These equations 
include, in addition to the inviscid terms, the leading viscous terms 
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according to Dunn's order-of -magnitude analysis. There are no dis- 
sipation terms in the simplified energy equation. 

The asymptotic method is supposed to solve the simplified equations 
with an error no larger than the error involved in dropping the missing 
viscous terms. It is obvious from Fig. 17*1 that the equations are 
better than the method used to solve them. The fact that the Lees- 
Reshotko neutral-stability curve is in goon agreement with the Dunn-Lin 
neutral curve also indicates that the main source of the failure of the 
asymptotic theory is not to be found in the equations used. Calculations 
at lover Mach numbers give good agreement of the three neutral curves at 

- 1.3, and fair agreement at = 1,6. It can be concluded that the 
asymptotic method is adequate for two-dimensional disturbances up tc 

- 1.6, and that the refinement of including the temperature fluctuations 
adds nothing to the adequacy of the theory. 

The experimental points shown in Fig. 17*1 were measured by Laufer 
and Vrebalovich (i 960 ). They are not important in the present context 
and will be discussed in Section 20.1. 

17.1.2 Effect of Mach number . 

The evolution of the neutral -stability curve with increasing Mach 
number is of considerable interest. Figure 17.2 gives neutral curves at 
five Mach numbers from 1.6 to 3 . 8 , The wave number is plotted against 
l/R to emphasize the high Reynolds number region. At = 1.6, the neutral 
curve still has the same character as in incompressible flow. There is 
a weak inviscid instability- and « approaches the inviscid neutral sub- 
sonic wave number c< st from above as R , In the sense of the 
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INVERSE REYNOLDS NUMBER, ^XIO 3 

Fig. 17.2 Effect of free-st ream Mach number on neutral-stability curve of 

insulated-wall boundary layer. 











definition given in the footnote on p. 75 > viscosity has a destabilizing 
effect with this type of neutral curve. 

At = 2.2, a local minimum appears in the upper branch of the 
neutral curve. This minimum indicates that the effect of viscosity has 
become stabilizing over a range of high Reynolds numbers. At - 2.6, 

this minimum is more pronounced and cx appears to approach from 

below. If this supposition is correct, then viscosity has a stabilizing 
effect from R up to the local minimum. At = 3,0, the minimum 

has virtually disappeared and oc clearly approaches o^, from below. 

The effect of viscosity is stabilizing, or almost so, over the entire 
range of Reynolds numbers. At the minimum has definitely 

disappeared and the wave number decreases monotonically along the. neutral 
curve from o( St to oi « 0. This neutral curve is of the same general type 
as the one in Fig. 9.2' (p, 7*0 far "the Falkner-Skan separation profile. 

The maximum instability is inviscid (at R -** oo), and the effect of; vis- 
cosity is only stabilizing. These statements can all be verified by 
studying the amplification rate as a function of Reynolds number, as will 
be done in Section 17.2.1. ■ • 

The neutral-stability curve at M-j - 4.2, plotted as e( vs JR, c h vs R 
and Cm ir* J lL t vs R, is shown in Fig. 17.3- At this Mach number there 
are two separate loops, and it is to be expected that the first, or lower, 
loop is the finite Reynolds number counterpart of the inviscid first mode 
and the second, or upper, loop is the counterpart of the inviscid second 
mode. The upper branch of the first loop tends to tx fiJ , as R « 

while the upper branch of the second loop tends to of 51 , c s . The lover 
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and on the 


branch of the first loop tends to o a, c 0 as 
upper loop it is not possible to assign a definite limit. These limits 
fit the identification of the two loops with the first and second modes. 

If the eigenfunctions of the neutral disturbances are examined, it is 
found that there is little phase change on the first loop, and a large 
phase change on the second loop. 

The shape of the neutral-stability curve of the second mode indicates 
that here too viscosity has only a stabilising influence. The much greater 
inviscid instability of the second mode as compared to the first mode is 


reflected in the smaller critical Reynolds number of the second mode. 

The lowest Mach number at which the second mode has been calculated 
(because of the c£ R limitation of the method of numerical integration) 
is = 3*6, where the critical Reynolds number is R - 830 , and the wave 
number is <?L = 0,306 and the phase velocity is cv = 0.904. 

Figure 17.4 shows the ot vs R neutral curves at - ^*5 and 4.8. 

At = 4.5 there are still two separate loops as at = 4.2, but at 
1*5^ = 4.8 an impoz'tant change has taken place. The two separate unstable 
regions have merged into a single unstable region enclosed by a single 
neutral curve. This single curve starts at the inviscid neutral sonic 
solution and ends at the inviscid neutral subsonic solution of the second 
mode. The factor which controls whether there are individual or separate 
neutral curves is whether <x ix is greater than or less than <* 4J . At 
* °<5i and there are two separate loops. At = 4.8, 

oC, x ot Sl and there is a single loop. 
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At every pi, £ In the interior of the single neutral curve there 
is an amplified solution. However, there are multiple viscous solutions, 
i.e., more than one c for a given , just as there are multiple 

Inviscid solutions, i.e., more than one c for a given oi , Hence at 

the same R where there is an amplified solution of one family, there 
can be a damped or neutral solution of another family. Some examples of 
multiple viscous solutions are shown in Figs. 19.4 and 19 . 5 . 

The neutral-stability curves of oi vs B and vs R at = 5.8 are 

given in Fig. 17.5 together with a few contours of = const. Although 
the merger of the two unstable regions is further advanced that 'at 

= 4.8, it is seen that these contours still form two distinct families. 
Consequently, it is still possible to identify disturbances as belonging 
to the first and second modes even after the merger. There is a- rapid 
change from first-mode to second -mode character at the apparent juncture 
of the two regions. The fact that both families of contours are - ' open to 
infinity is another indication that the effect -of viscosity is purely 
destabilizing. . 1 
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Fig. 17 • 5 Neutral-stability curvet of photo velocity and wove number and contour* of amplified 

disturbance* at M, = 5.8 
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Fig. 17.11 Amplitude ratio vs dimensioniess frequency for 
four Reynolds numbers at.Mx = 3.8 
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•^ig* 17.3-2 Frequency response of insulated-wall boundary layer at M “ 5.8 
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17*13 Amplitude ratio of most unstable disturbance at R = 1500 as a function 0 


free-stream Mach number. Insolated-vrall boundary layer. 
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Pig. 17*13 Amplitude ratio of most unstable disturbance at R - 1500 as a function of 

free-stream Mach number. Insulated-wall boundary layer. 
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Fig. 18.2 Effect of wave angle on the amplification 
rate of first-mode three-dimensional 
disturbances at R = 1500 
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R 


$ difference 


M-, 


6th order 8th order 

C* ^ «C; X |() 3 0< C; it |0* 


1.3 

500 

0.075 

45° 

0.883 

0.824 

7.2 

1.3 

1500 

O.OoO 

45° 

1.467 

1.445 

1.5 

1.6 

500 

0.070 

55° 

0.974 

0.874 

11.4 

1.6 

1500 

0.050 

55° 

1.384 

1.346 

2.8 

2.2 

500 

0.055 

60° 

1.198 

1.066 

12.4 

2.2 

800 

0. Oh-5 

60° 

1.391 

1.300 

7.0 

2.2 

1500 

0.035 

6o° 

1 . 32s 

1.273 

4.3. 

4.5 

500 

0.045 

6 o° 

1.117 

1.039 

7.5 

4.5 

1500 

0.050 

60° 

1.641 

1.613 

1.7 

5.8 

500 

0.050 

0 

55 

0.790 

0.736 

7.3 

5.8 

1500 

0.060 

55° 

1.403 

1.384 

1.4 

10.0 

1500 

0.040 

55° 

0.444 

0.434 

2.3 


Table 18.1 Comparison of amplification rates for three-dimensional 
disturbances as computed from sixth-order and eighth- 
order systems of equations at several Mach numbers. 




$ change in « c ; 


zero term 

M 1 *= 2.2 


10.0 

x, £* - 
*1 ’7 

+ 1.4 

+ 0.7 

+ 0.9 

2 Jfc »v' 

rt 

0 

- 0.1 

- 0.4 

k.( iU \ T ' 

- 0.1 

- 1.2 

- 2.4 

^ 3 J1 +** ^ 

+ 4.1 

+ 1.4 

+ .2.3 

All zero 

+ 5.3 

+ 0.8 

+ 0.4 


Table 18.2 Comparison of individual energy-equation 

dissipation terms on amplification rate at 
R = 1500 for three Mach numbers. 
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Fig. 19.2 Neufral-stabiliiy curve of wave number a? . 

A/J ! = 5.8 (cooled boundary layer; 6 W = 0.10} 
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Fig. 19.3 Neutral-stability curve of wave number at 
Mi = 5.3 (cooled boundary layer; 0 W = — 0.1 0) 
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Fig. 19.6 Phase velocity vs wave number for cooled-walf 
boundary layer at Mi = 5.0, t w /T r = 0/25 





fi e . 19.7 EFFECT OF COOLING ON MAXIMUM AMPLIFICATION RATES M| =5.8 
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Fig. 19.8 Amplification rave or the most unstable 
frequency at R = 1500 vs Reynolds Number for 
the insulated and cooled boundary Sayers 
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Fig. 19*9 Effect of cooling at M = 5. S on the amplification at R = 1500 of two 

types of disturbances. 
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